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CHAPTER  1:  INTRODUCTION 

Fingerprint  identification  is  one  of  the  principal  methods 
used  by  criminal  investigators  to  obtain  evidence.   The 
uniqueness  of  individual  fingerprints  allows  for  positive 
identification  of  prints  left  at  the  scene  of  a  crime. 
Prints  cannot  be  altered;  thus,  they  provide  a  permanent 
form  of  identification. 

Fingerprints  can  be  obtained  from  different  types  of  sur- 
faces.  One  type  of  surface  which  presents  difficulty  in 
obtaining  prints  is  human  skin.   The  difficulty  arises  from 
physiological  processes  which  occur  in  the  body  that  result 
in  rapid  deterioration  of  prints  on  skin.   Therefore,  prints 
on  human  skin  must  be  recovered  quickly  before  they  become 
unrecognizable.   In  a  situation  where  a  victim  has  been 
traumatized,  such  as  rape,  a  period  of  time  usually  passes 
before  the  crime  is  reported.   As  the  time  following  place- 
ment of  the  print  lengthens,  it  becomes  increasingly  diffi- 
cult for  investigators  to  obtain  distinguishable  prints. 

The  primary  purpose  of  this  project  is  to  increase  the 
recovery  time  for  prints  obtained  from  human  skin  through 
digital  image  processing  methods.   Initially,  pictures  of 
prints  are  recorded  as  digital  images.   The  image  is  then 
processed  using  image  enhancement  techniques.   Through  these 


techniques,  it  is  possible  to  enhance  the  ridge  structure  of 
degraded  fingerprint  images  in  order  to  produce  a  recogniz- 
able print. 

There  is  a  wide  variety  of  image  enhancement  techniques 
available  for  improving  image  quality.   The  features  to  be 
enhanced  determine  which  type  of  techniques  are  best  suited 
for  the  task.   Most  routines  attempt  to  accentuate  areas  of 
differing  contrast.   For  instance,  highlighting  edges 
between  contrasting  regions  is  one  method  often  used.   Many 
times  a  combination  of  several  techniques  will  produce  the 
best  results.   Examination  of  the  statistics  for  an  image 
can  often  provide  insight  as  to  the  effectiveness  of  the 
enhancement  methods. 


CHAPTER  2:  FINGERPRINT  IDENTIFICATION  AND  RECOVERY 

The  term  dermatoglyphics  was  coined  around  1926  to  describe 
the  science  of  the  study  of  skin  patterns  [3].  It  particu- 
larly applies  to  the  patterns  of  the  specialized  skin  of  the 
inferior  surfaces  of  the  hands  and  feet,  those  being  the 
soles  of  the  feet  and  the  palms  and  fingertips  of  the  hand. 
Dermatoglyphics  is  concerned  primarily  with  the  morphologi- 
cal and  physiological  aspects  of  these  skin  patterns. 

The  skin  of  the  inferior  surfaces  of  the  hands  and  feet  is 
unlike  that  of  any  other  part  of  the  body.   These  areas, 
which  contain  neither  hair  nor  sebaceous  (oil)  glands,  pos- 
sess an  abundance  of  relatively  large  sweat  glands.   Perhaps 
the  most  important  characteristic  of  the  skin  in  these  areas 
is  the  continuously  corrugated  surface  with  narrow  ridges. 

One  important  practical  application  of  dermatoglyphics  is  in 
the  area  of  anthropometric  identification,  specifically, 
fingerprint  identification.   A  fingerprint  is  formally 
defined  as  a  reproduction  of  the  pattern  or  design  formed  by 
the  ridges  of  the  first  joint  of  a  finger  or  thumb.   To 
date,  no  two  fingerprints  have  been  found  to  be  exactly 
alike.   Fingertip  patterns  are  formed  during  the  early 
embryonic  stage  and  never  change  during  the  entire  course  of 
a  human  life.   Even  injury  does  not  change  the  patterns 


since  the  ridges  are  known  to  heal  exactly  as  before  to 
within  a  millionth  of  an  inch.   Many  torturous  attempts  to 
alter  ridge  patterns  have  failed  due  to  the  natural  healing 
process. 

The  pattern  arrangements  of  fingertip  ridges  have  been 
noticed  by  man  for  centuries.   It  is  known  that  the  Chinese 
used  fingerprints  to  sign  documents  dating  back  to  the  third 
century  B.C.   A  Chinese  contract  of  loan  executed  nearly 
1200  years  ago  contains  the  imprint  of  a  fingertip  in  a  pat 
of  clay  which  is  attached  to  the  document.   An  inscription 
on  the  document  states,  "The  two  parties  have  found  this 
just  and  clear  and  have  affixed  impressions  of  their  fingers 
to  serve  as  a  mark." 

Scientific  descriptions  of  the  patternings  of  fingers  and 
palms  date  back  to  the  late  seventeenth  century.   The  first 
item  in  modern  literature  relating  to  fingerprint  identifi- 
cation was  published  by  Henry  Faulds  in  1880  in  an  English 
journal  called  Nature.   In  this  short  article  Faulds  points 
out  that  chance  fingerprints  left  at  the  scene  of  a  crime 
would  provide  positive  identification  of  the  offender  when 
apprehended. 

The  modern  system  of  fingerprint  identification  originated 
in  the  late  nineteenth  century  from  work  done  by  Faulds,  Sir 
Francis  Galton,  Sir  Edward  Richard  Henry,  and  many  others. 


The  current  system  is  called  the  Henry  system  after  Sir 
Henry.   In  the  Henry  system,  fingerprint  characteristics  are 
divided  into  three  general  patterns  called  loops,  arches, 
and  whorls  [2].   Prints  are  categorized  by  the  relative 
locations  of  specific  patterns.   The  patterns  are  described 
by  the  number  of  ridges  which  comprise  them.   An  example  of 
a  whorl  is  shown  in  Figure  2-1. 


Bifurcation 


Ridge  Ending 


Figure  2-1  Example  of  a  whorl  with  bifurcation 
and  ridge  ending 

The  most  important  features  of  fingerprints  are  bifurca- 
tions and  ridge  endings  [8].   Bifurcations  are  locations  in 
a  fingerprint  where  ridges  branch  out  in  different  direc- 
tions.  Ridge  endings  are  simply  locations  where  a  specific 
ridge  terminates.   Identification  of  a  print  is  accomplished 
by  pinpointing  the  exact  locations  of  points  corresponding 
to  bifurcations  and  ridge  endings.   Each  state  has  its  own 
laws  governing  the  exact  number  of  points  required  to  make  a 
positive  identification.   For  example,  the  State  of  Kansas 


requires  nine  such  points  whereas  California  requires  six- 
teen points  to  make  a  positive  identification  of  a  print. 
Defects  such  as  permanent  scars  are  also  used  in  fingerprint 
identification  since  they  render  a  print  unique. 

The  most  difficult  task  confronting  criminal  investigators 
is  the  recovery  of  fingerprints  from  the  scene  of  a  crime. 
Fingerprints  left  at  the  scene  of  a  crime  are  basically  of 
three  types:  visible,  plastic,  and  latent  prints.   Visible 
prints  are  transferred  from  smearings  of  the  hand  with  ink, 
paint,  blood,  or  other  substances  which  come  into  contact 
with  a  surface  of  contrasting  color.   Plastic  prints  are 
formed  from  ridge  impressions  in  pliable  materials  such  as 
soap,  putty,  wax,  and  dust.   The  most  difficult  type  of 
prints  to  recover  are  latent  prints.   These  are  ridge 
impressions  which  are  formed  on  surfaces  by  perspiration  or 
oil  from  the  fingers.   Latent  prints  are  scarcely  visible  so 
development  of  the  print  is  necessary  for  further  study. 

Many  methods  are  used  to  recover  latent  fingerprints.   The 
most  common  method  invloves  dusting  a  surface  with  a  powder 
of  contrasting  color.   This  method  works  well  on  hard,  non- 
porous  surfaces.   Another  popular  method  is  the  treatment  of 
a  surface,  such  as  paper,  with  a  silver  nitrate  solution. 
There  are  many  other  esoteric  methods  such  as  interferomet- 
ric  and  infrared  microscopy,  laser  illumination,  metal 


evaporation/deposition,  and  x-ray  electronography. 

Human  skin  is  one  surface  from  which  it  is  extremely  diffi- 
cult to  recover  latent  prints.   A  print  left  on  the  surface 
of  human  skin  will  deteriorate  after  a  period  of  time  due  to 
perspiration  and  oil  secretion  from  beneath  the  skin.   Thus, 
the  recovery  time  for  a  latent  print  left  on  human  skin  is 
limited.   The  actual  recovery  time  is  dependent  upon  the 
persons  involved  because  oil  and  skin  compositions  differ. 
Good  prints  have  been  recovered  from  live  persons  within  a 
period  of  ninety  minutes  and  from  cadavers  up  to  forty-one 
hours  after  placement. 

The  recovery  of  latent  prints  from  human  skin  is  a  great 
asset  to  criminal  investigators  involved  in  cases  such  as 
strangulation  and  rape.   The  most  widely  used  method  is  the 
iodine-silver  plate  method.   It  has  been  adopted  by  the 
F.B.I,  as  well  as  many  State  and  Federal  agencies.   The 
prints  for  this  project  were  obtained  using  a  modification 
of  this  method  [5]. 

The  modified  iodine-silver  plate  method  begins  with  the 
preparation  of  the  recovery  medium  which  is  a  two-inch  piece 
of  35  mm  camera  film.   This  is  immersed  in  a  fixer  solution 
to  strip  away  the  silver  salts  so  that  only  transparent  ace- 
tate remains.   The  clear  acetate  is  dipped  in  a  solution  of 
leuco  crystal  violet  in  xylene.   Leuco  crystal  violet  is  a 


reagent  used  to  detect  iodine. 

Once  the  treated  acetate  is  dry,  the  skin  area  is  prepared 
for  the  print  transfer.   The  suspected  area  is  fumed  with 
iodine  vapor.   Iodine  vapor  is  known  to  be  absorbed  by  the 
body  oils  that  generally  comprise  the  ridges  of  latent 
prints.   Continued  fuming  of  the  area  reveals  the  ridge 
detail  of  the  print.   When  the  print  is  visible,  the  treated 
acetate  is  pressed  over  the  print  with  uniform  pressure 
until  transferred.   A  thin  coating  of  rubber  cement  is  then 
brushed  over  the  acetate  to  preserve  the  print. 


CHAPTER  3:  ENHANCEMENT  TECHNIQUES 

3.1  IMAGE  DEFINITIONS 

Image  enhancement  routines  process  digital  images  in  order 
to  refine  certain  features.   A  digital  image  may  be  repre- 
sented as  a  two-dimensional  light-intensity  function  denoted 
as  f(x,y)  where  x  and  y  are  independent  variables  denoting 
spatial  coordinates.   Specific  coordinates  are  referred  to 
as  picture  elements  or  pixels.   The  value  of  f  describes  the 
brightness  of  the  image  at  the  coordinates  (x,y).   This 
brightness  intensity  is  divided  into  discrete  values  called 
gray  levels.   A  digital  image  may  be  considered  as  a  matrix. 
Each  element  of  the  matrix  represents  the  gray  level  value 
of  the  pixel  at  that  location.   The  location  of  a  pixel  in 
the  matrix  is  identified  by  its  spatial  coordinates. 

Enhancement  operations  may  be  represented  as  mathematical 
transformations.   In  general,  the  digital  image  f(x,y)  is 
transformed  to  an  image  g(x,y)  according  to  a  specific  math- 
ematical rule.   A  variety  of  rules  may  be  derived  to  enhance 
specific  features  contained  in  an  image. 

Enhancement  techniques  are  generally  classified  as  either 
spatial-domain  or  frequency-domain  methods  [6].   Spatial- 
domain  methods  directly  manipulate  the  pixels  of  an  image. 
Frequency-domain  methods  modify  spatial  frequency  informa- 
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tion  to  enhance  corresponding  features  in  the  spatial 
domain.   The  Fourier  transform  is  used  to  transform  the  ori- 
ginal image  to  the  freguency  domain. 

3.2   SPATIAL-DOMAIN  METHODS 

Spatial-domain  methods  work  on  a  variety  of  principles.   The 
most  elementary  operation  is  known  as  thresholding.   In  this 
method,  a  cutoff  value  called  a  threshold  is  specified. 
Gray  levels  either  above  or  below  this  threshold  are  changed 
to  a  specified  value  by  the  transformations 


g(x,y)  = 


or 


g(x,y)  = 


f(x,y)  if  f(x,y)  >  t 

c  f(x,y)  <  t 

f(x,y)  if  f(x,y)  <  t 

c  f(x,y)  >  t 


(3.2-1) 


where  c  is  a  specified  gray  level  and  t  is  the  threshold 
value.   This  method  is  effective  when  a  large  degree  of  con- 
trast is  evident  in  the  original  image. 

Thresholding  is  actually  a  specific  case  of  a  technigue 
known  as  contrast  stretching.   This  is  a  technigue  where 
gray  levels  from  f(x,y)  are  mapped  to  g(x,y)  according  to  a 
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specific  transformation  function.   Let  r  and  s  represent 
gray  level  values  in  f(x,y)  and  g(x,y),  respectively.   The 
transformation  function  is  given  by  s  =  T(r)  where  T  is  a 
mathematical  operator  defined  on  r.   A  transformation  func- 
tion is  constrained  to  be  a  continuous,  monotonically 
increasing  function  whose  range  is  equivalent  to  its  domain. 
Examples  of  transformation  functions  are  shown  graphically 
in  Figure  3.2-1. 


255 


255 


255 


Figure  3.2-1  Gray-level  transformation  functions 

Another  class  of  spatial-domain  methods  are  known  as  histo- 
gram modification  techniques  [6].   A  histogram  is  a  plot 
which  provides  a  description  of  the  distribution  of  gray 
levels  in  an  image.   Gray-level  distributions  are  described 
mathematically  using  a  probability  density  function.   The 
probability  density  function  is  defined  according  to  the 
relative  frequency  of  gray  levels  appearing  in  an  image 
expressed  as 
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Pr(r)  =  t  /  N  (3.2-2) 

where  t  is  the  number  of  pixels  with  gray  level  r  and  N  is 
the  total  number  of  pixels  in  f(x,y). 

Histogram  modification  techniques  map  values  by  specifying 

the  probability  density  function  for  the  enhanced  image. 

Let  p  (s)  denote  a  specific  probability  density  function  for 
s 

the  enhanced  image  g(x,y)  and  s  =  T(r)  denote  the  transfor- 
mation function.   Histogram  modification  techniques  are 
based  on  the  relationship 

Ps(s)  =  [pr(r)dr/ds]|r=T"1(s)         (3.2-3) 

where  T   (s)  is  the  inverse  transformation  function. 

The  most  common  histogram  modification  technique  is  known  as 
histogram  equalization.  Enhancement  is  accomplished  through 
equalization  using  the  transformation  function 

s  =  T(r)  =   |  p  (w)dw  (3.2-4) 


'  Pr(wl 


0 

where  w  is  a  dummy  variable.   This  transformation  attempts 
to  distribute  the  gray  levels  from  the  density  function 
p  (r)  as  a  uniform  density  function,  denoted  as  p  (s),  in 
order  to  effect  an  increase  in  the  dynamic  range  of  the 
pixels  in  the  image. 


12 


A  useful  variety  of  enhancement  operations  are  image  sharp- 
ening techniques  [6].   These  are  employed  primarily  as  a 
means  for  highlighting  edges  in  an  image.   In  the  spatial 
domain,  differentiaton  is  the  operation  used  to  emphasize 
edge  detail.   An  edge  is  a  boundary  between  contrasting 
regions  characterized  by  an  abrupt  change  in  the  gray  levels 
of  adjacent  pixels. 

Differentiation  is  usually  executed  using  the  gradient  oper- 
ator.  The  gradient  of  f(x,y)  is  a  vector  defined  as 


G[f(x,y)]  = 


^x 


(3.2-5) 


The  gradient  is  a  directional  derivative  which  points  in  the 
direction  of  the  maximum  rate  of  increase  of  f(x,y). 

Gradients  are  implemented  through  the  use  of  spatial  masks. 
A  mask  is  a  small  two-dimensional  array  with  coefficients 
which  are  chosen  to  recognize  a  certain  property  in  an 
image.   A  spatial  mask  is  a  local  operator  in  that  it  is 
only  concerned  with  a  single  point  and  its  surrounding 
neighborhood.   A  3  x  3  spatial  mask  with  corresponding  pixel 
locations  is  shown  in  Figure  3.2-2. 
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f(x-l,y-l) 

w2 
f(x-l,y) 

w3 
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f(x,y+l) 

f(x+l,y-l) 

w8 
f(x+l,y) 

w9 

f(x+l,y+l) 

Figure  3.2-2   3x3  spatial  mask 


Consider  the  3x3  mask  and  image  neighborhood  as  row- 
ordered  column  vectors  w  and  f,  respectively.   These  vectors 
are  written  as 


f  = 


- — 

f(x-l, 

y-1) 

f(x-l, 

t 

y+l) 

. 

f(x+l, 

y-D 

f(x+l'( 

y+l) 

w  = 


"~~ 

wl 

• 

. 

_"3_ 

• 

_L_ 

^7 

* 

« 

W9 

i —      1 

(3.2-6) 


The  computation  of  the  mask  operation  is  executed  as  the 
inner  product  of  the  two  vectors: 

w'f  =  Wjf(x-l,y-l)  +  w2f(x-l,y)  +     +  wgf(x+l,y+l) 

(3.2-7) 
The  value  f(x,y)  is  replaced  by  the  magnitude  of  the  inner 
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product  to  produce  the  gradient. 

If  appropriate  masks  are  used,  the  magnitude  of  the  gradient 
reveals  points  whose  neighbors  have  greatly  differing  gray 
levels.   These  points  correspond  to  singularities  in  f(x,y). 
Edges  are  represented  as  a  connected  boundary  of  singulari- 
ties.  Many  types  of  masks  have  been  derived  which  locate 
edges  oriented  in  various  directions  [9]. 

Image  smoothing  operations  are  used  to  remove  noise  which 
may  be  present  in  an  image.   One  such  method  is  known  as 
neighborhood  averaging.   This  is  performed  by  replacing  a 
pixel  with  the  average  of  the  surrounding  neighborhood  using 
a  spatial  mask.   A  mask  used  for  averaging  over  a  3  x  3 
neighborhood  is  shown  in  Figure  3.2-3. 
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Figure  3.2-3  Spatial  mask  used  for 
neighborhood  averaging 


Median  filtering  is  another  smoothing  technique.   A  median 
filter  replaces  a  pixel  with  the  median  of  the  surrounding 
neighborhood.   The  advantage  of  this  method  over  neighbor- 
hood averaging  is  a  reduction  in  blurring  which  preserves 
edge  content  in  an  image  [ 7 ] . 
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Psuedo-color  processing  is  a  relatively  new  class  of  spa- 
tial-domain methods.   Many  operations  are  possible  within 
this  powerful  area  of  enhancement  techniques.   Psuedo-color 
techniques  are  based  on  the  principle  that  the  human  eye  has 
the  ability  to  discern  a  wide  range  of  colors  and  intensi- 
ties.  This  ability  is  limited  when  viewing  monochrome 
images . 

The  simplest  operation  in  psuedo-color  techniques  is  called 
density  slicing.   Density  slicing  is  useful  in  displaying 
the  spatial  locations  of  specific  gray  levels.   The  basic 
idea  behind  this  process  is  to  assign  a  color  to  each  pixel 
based  on  its  intensity.   This  is  analogous  to  slicing  the 
two-dimensional  image  function  f(x,y)  with  a  plane  and 
assigning  a  color  to  that  plane.   A  geometric  interpretation 
of  this  process  is  shown  in  Figure  3.2-4. 
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f(x,y)  A 


Figure  3.2-4  Geometrical  interpretation  of  density  slicing 


17 


3.3  FREQUENCY-DOMAIN  METHODS 

Frequency-domain  methods  enhance  images  by  modifying  spatial 
frequencies  which  correspond  to  particular  features  in  the 
spatial  domain.   Generally,  images  are  transformed  from  the 
spatial  domain  to  the  frequency  domain  via  the  discrete 
Fourier  transform  (DFT). 

For  an  M  x  N  image,  the  two-dimensional  DFT  is  computed  as 

M-l  N-l 
F(u,v)  =  (1/MN)   V~  Y~  f  (x,y)exp[-j2Tt(ux/M  +  vy/N)  ] 

x=0  y=0 

(3.3-1) 
for  the  spatial  frequencies  u  =  0,...,M-1  and  v  =  0,...,N-1. 
The  DFT  is  a  separable  transform,  so  the  two-dimensional  DFT 
can  be  computed  using  a  one-dimensional  FFT  algorithm  [6]. 
This  reduces  the  number  of  operations,  which  in  turn  pro- 
vides a  significant  increase  in  computing  speed  [1]. 

Filtering  is  the  basic  operation  used  in  frequency-domain 
methods.   Frequency  information  is  characterized  by  its  dis- 
tance from  the  origin  of  the  frequency  plane.   High  spatial 
frequencies  are  located  furthest  away  from  the  origin  while 
low  frequencies  are  closest  to  the  origin.   Filtering  mod- 
ifies spatial  frequencies  according  to  their  distance  from 
the  origin  in  the  frequency  plane.   The  extent  to  which  the 
frequency  values  are  modified  is  defined  by  a  transfer  func- 
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tion.   Values  of  this  function  are  restricted  to  the  inter- 
val [0,1].   A  cutoff  point  corresponding  to  a  radius  about 
the  origin  determines  the  rate  of  transition  from  the  min- 
imum to  the  maximum  values  of  the  transfer  function.   The 
value  of  the  transfer  function  at  the  cutoff  point  is  gener- 
ally equal  to  half  of  its  maximum  value. 

Let  G(u,v)  denote  the  function  for  the  enhanced  image  in  the 
frequency  domain.   Consider  a  transfer  function  in  the  freq- 
uency domain  given  by  H(u,v).   The  filtering  operation  is 
computed  as 

G(u,v)  =  F(u,v)H(u,v)  (3.3-2) 

This  operation  corresponds  to  convolution  in  the  spatial 
domain.   Transforming  the  filtered  image  back  to  the  spatial 
domain  yields  the  enhanced  image  g(x,y). 

Image  sharpening  is  accomplished  in  the  frequency  domain 
through  the  implementation  of  highpass  filters.   High-freq- 
uency components  in  the  frequency  domain  correspond  to  edge 
content  in  the  spatial  domain.   Edges  are  accentuated  by 
attenuating  the  low-frequency  components  using  a  highpass 
transfer  function.   High-frequency  components  may  be  further 
emphasized  by  adding  a  constant  to  the  highpass  transfer 
function.  That  is, 

H  =  H  +  c  (3.3-3) 

where  c  is  a  constant. 
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CHAPTER  4:  PROJECT  DESCRIPTION 

4.1  IMAGE  PROCESSING  SYSTEM 

The  objective  of  this  project  is  to  enhance  degraded  finger- 
print images  using  digital  processing  methods.   This  re- 
guires  that  pictures  of  the  prints  be  recorded,  via  a  digi- 
tizing camera,  as  digital  images  .   Subseguently,  the  images 
may  be  viewed  on  a  video  display  monitor.   These  operations 
are  administered  via  a  general-purpose  computer.   The  cru- 
cial link  among  all  of  the  aforementioned  components  is  the 
image  processing  system.   An  image  processing  system  con- 
tains storage  modules  called  frame  buffers  which  store  the 
images  during  the  acguisition  and  display  stages. 

A  Grinnell  GMR  270-series  image  processing  system  was  used 
to  digitize  and  display  pictures  of  prints.   These  pictures 
are  captured  as  monochrome  images,  using  a  television  camera 
as  the  input  device.   Images  are  generated  as  a  dot  matrix 
grid  consisting  of  512  lines  with  512  picture  elements 
spaced  along  each  line.   The  Grinnell  provides  a  dynamic 
brightness  range  of  256  discrete  gray  levels  contained  in 
the  interval  [0,255]  for  each  pixel.   A  value  of  0  corre- 
sponds to  the  dark  elements  while  255  corresponds  to  the 
brightest  elements.   Images  are  displayed  on  a  Mitsubishi 
Model  #C3922  LPK  high-resolution  color-television  monitor. 
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A  DEC  VAX  11/750  digital  computer  was  used  as  the  host  com- 
puter to  control  operation  of  the  Grinnell. 

4.2  COMPUTER  PROGRAMS 

Two  programs  were  written  to  process  the  degraded  print 
images:  ENHANCE  and  COLOR.   Both  programs  are  interactive, 
which  provides  the  user  with  some  degree  of  flexibility. 
They  are  also  menu-driven,  facilitating  ease  of  operation. 

ENHANCE  is  an  interactive  program  which  performs  a  variety 
of  spatial-domain  and  frequency-domain  enhancement  routines. 
A  collection  of  unary  and  binary  arithmetic  operations  are 
included  so  that  various  mathematical  manipulations  are  pos- 
sible.  Statistical  operations  are  incorporated  to  allow 
examination  of  the  distribution  of  gray  level  values  in  an 
image. 

The  user  is  required  to  provide  values  of  certain  parameters 
in  order  to  tailor  the  routines  to  specific  applications. 
Menus  are  used  to  prompt  the  user  for  these  values.   Each 
menu  displays  the  options,  current  values  of  parameters,  and 
the  last  operation  performed.   The  program  records  up  to 
fifteen  operations  in  a  memory  stack  to  assist  the  user  in 
recalling  previous  operations. 

ENHANCE  contains  two  working  registers,  X  and  Y,  and  two 
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storage  registers,  1  and  2.   Each  register  is  a  512  x  512 
complex- valued  data  array.   Unary  operations  are  performed 
on  data  in  the  X  register.   Binary  operations  reguire  data 
in  both  the  X  and  Y  registers.   Results  are  always  returned 
to  the  X  register.   Data  which  was  previously  stored  in  the 
X  register  is  not  saved  due  to  the  cost  in  speed  and  memory. 
The  storage  registers  are  available  for  storing  intermediate 
results. 

External  data  files  are  input  to  the  program  as  real  values 
between  0  and  255.   The  imaginary  component  is  set  to  zero 
upon  input.   Data  files  are  written  as  integer  values  be- 
tween 0  and  255.   Only  the  real  component  of  the  X  register 
is  written  to  external  data  files. 

COLOR  is  a  program  used  to  perform  psuedo-color  enhancement, 
specifically,  density  slicing.   This  is  a  technigue  in  which 
gray  levels  are  divided  into  intervals  and  a  specific  color 
assigned  to  each  interval.   The  image  is  then  displayed 
using  color  planes  to  show  the  spatial  locations  of  specific 
gray  levels. 

COLOR  contains  four  registers.  The  registers  are  512  x  512 
integer-valued  data  arrays.  The  main  register  contains  the 
image  to  be  operated  upon.  Red,  green,  and  blue  color  reg- 
isters are  used  to  display  the  partitioned  intervals.  Data 
may  be  accessed  from  any  of  the  four  registers. 
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CHAPTER  5:  RESULTS 

5.1  INTRODUCTION 

This  project  involves  the  enhancement  of  degraded  finger- 
print images.   Ridge  detail  is  the  most  important  feature  in 
fingerprints  used  for  identification.   Therefore,  the  main 
thrust  of  the  enhancement  technigues  used  was  to  enhance  the 
contrast  between  the  regions  corresponding  to  the  ridges  and 
valleys  of  the  fingerprints. 

Samples  for  this  project  were  provided  by  Jeff  Payne  of  the 
Department  of  Chemistry  at  Kansas  State  University.   The 
prints  appeared  as  a  faint  bluish  outline  on  a  piece  of 
transparent  35  mm  film  acetate.   Samples  were  the  results  of 
prints  placed  on  a  subject's  arm  and  lifted  at  regular  in- 
tervals. 

A  television  camera  connected  to  the  Grinnell  image  process- 
ing system  was  used  to  capture  the  samples  as  monochrome 
digital  images.   The  samples  were  placed  on  a  light  table  in 
order  to  highlight  the  prints  while  photographing.   Programs 
on  the  VAX  11/750  contained  in  the  KSU  IMAGE  library  were 
utilized  to  control  operation  of  the  television  camera  and 
frame  buffers.   The  keyword  IMAGE  must  be  entered  to  ini- 
tialize the  programs.   Once  initialized,  the  keyword  CAMERA 
is  used  to  activate  the  television  camera.   When  the  desired 
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image  is  displayed  on  the  monitor,  the  keyword  STOIMG  trans- 
fers the  image  to  a  512  x  512  data  file.  Data  files  require 
approximately  514  blocks  of  memory;  thus,  care  must  be  exer- 
cised so  that  the  user's  disk  quota  is  not  exceeded. 

The  sample  documented  in  this  thesis  was  recovered  two  hours 
after  placement.   Figure  5.1-1  shows  a  picture  of  this  sam- 
ple.  A  visual  analysis  of  this  image  indicates  some  faint 
ridge  detail  in  the  center  of  the  image,  forming  a  pattern 
which  appears  to  be  a  loop.   Degradation  is  most  severe  in 
the  lower  portion  of  the  image  to  the  extent  that  individual 
ridges  are  indistinguishable.   Ridges  in  the  left  half  of 
the  image  are  extremely  faint.   The  bubbles  which  appear  in 
the  image  are  a  result  of  the  recovery  process  and  cannot  be 
removed . 

The  program  COLOR  was  used  to  implement  density  slicing 
techniques  as  an  aid  in  determining  which  gray  levels  com- 
prise the  ridge  detail  in  the  image.   This  was  done  by 
dividing  the  gray  scale  into  intervals  and  assigning  a  color 
to  each  partitioned  interval.   A  display  of  the  composite 
image  reveals  the  interval  to  which  a  spatial  coordinate 
belongs. 

Ridges  were  found  to  lie  in  the  gray- level  range  from  6  to 
160,  with  the  majority  of  the  distinguishable  ridges  con- 
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tained  in  the  range  from  100  to  120.   Ridge  information  in 
the  lower  section  of  the  image,  where  the  worst  degradation 
occurs,  falls  into  the  interval  from  6  to  100.   The  far  left 
and  right  regions  of  the  image  reveal  very  faint  ridges 
which  occupy  the  range  from  130  to  160.   The  one  character- 
istic which  became  evident  through  density  slicing  was  that 
the  ridges  were  not  uniform  with  many  different  gray  values 
comprising  a  single  ridge. 

A  statistical  analysis  of  the  gray  level  content  in  the 
image  was  performed  using  the  program  ENHANCE.   Gray  levels 
for  this  image  were  found  to  lie  in  the  range  from  6  to  252. 
These  gray  levels  are  distributed  according  to  the  histogram 
shown  in  Figure  5.1-2.   A  mean  value  of  135  with  a  variance 
of  758  was  calculated  from  this  distribution.   The  histogram 
indicates  that  most  of  the  gray  level  information  is  con- 
tained in  the  interval  from  60  to  245. 

Enhancement  operations  were  executed  through  the  program 
ENHANCE.   Various  combinations  of  spatial-domain  and 
freguency-domain  operations  were  tried  until  the  best  tech- 
nigues  were  discovered.   No  single  technigue  is  optimal 
since  the  range  of  contrast  and  brightness  can  vary  greatly 
within  any  given  image.   Filtering  in  the  freguency  domain 
was  found  to  produce  the  best  results  for  sharpening  ridge 
detail.   Spatial  domain  technigues  proved  useful  in  enhanc- 
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ing  contrast  in  the  image.   Two  procedures  which  yielded 
good  results  are  detailed  in  this  chapter. 

5.2  PROCEDURE  ONE 

The  first  procedure  utilizes  contrast  stretching  operations 
to  increase  the  degree  of  contrast  present  in  the  image. 
Contrast  stretching  functions  well  when  most  of  the  gray 
levels  are  contained  within  a  well-defined  range.   Filtering 
is  also  used  to  enhance  the  edges  between  the  ridges  and  the 
valleys. 

A  linear  transformation  function  was  used  initially  so  that 
ridge  detail  would  be  preserved.   The  transformation  is 
given  as 

0   if       f(x,y)  <  c. 


g(x,y)  =  < 


255[f(x,y)  -  c±]  c1   <   f(x,y)  <  ch 

(ch  -  Cl) 

255  f(x,y)  >  ch 


(5.2-1) 
where  c,   and  c,   are  the  cutoff  points.   This  transforma- 
tion is  shown  graphically  in  Figure  5.2-1.   The  cutoff 
points  c,  and  c,  chosen  for  the  image  in  Figure  5.1-1  were 
70  and  240,  respectively.   These  cutoff  points  caused  a 
small  percentage  of  pixels  outside  of  the  cutoff  range  to  be 
transformed  to  the  extreme  values  0  and  255.   Figure  5.2-2 
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shows  a  picture  of  the  transformed  image.   Although  the 
improvement  in  contrast  is  not  glaringly  obvious  to  the 
human  eye,  it  is  guite  beneficial  to  the  subseguent  filter- 
ing operations.   A  plot  of  the  histogram,  shown  in  Figure 
5.2-3,  exhibits  an  expanded  gray  level  range  which  gives 
credence  to  the  name  contrast  stretching.   The  increase  in 
the  degree  of  contrast  is  reinforced  by  the  increase  in  var- 
iance to  1904.   A  reduction  in  the  mean  value  to  103 
occurred  as  a  result  of  the  transformation. 


f(x,y) 


Figure  5.2-1  Linear  contrast  stretching 
transformation  function 


Filtering  is  performed  to  accentuate  the  edges  between  the 
ridges  and  valleys  of  the  print.   This  operation  requires 
that  the  origin  in  the  frequency  domain  be  translated  to  the 
center  of  the  frequency  plane.   Translation  is  achieved 
through  the  relationship 

f  (x,y)exp[j2ir(u0x  +  vQy)  ]  <=>  F(u-u0,v-vQ) 

(5.2-2) 
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where  the  symbol  <=>  denotes  "Fourier  transform  pair." 
Once  translated,  the  image  is  transformed  to  the  frequency 
domain  via  the  DFT.   A  useful  relationship  for  displaying 
the  transformed  image  is  given  by 

Ff (u,v)  =  Log[l  +  |  F(u,v)  |]  (5.2-3) 

This  relationship  is  useful  for  displaying  pixels  of  rela- 
tively small  magnitude.   A  picture  of  F'(u,v)  is  shown 
in  Figure  5.2-4.   Examination  of  this  picture  indicates  that 
most  of  the  frequency  information  is  located  towards  the 
center  of  the  frequency  plane.   A  significant  amount  of 
information  lies  on  the  horizontal  and  vertical  axes.   The 
spikes  located  on  the  diagonals  correspond  to  noise  and 
periodicity  resulting  from  circular  symmetry  in  the  ridge 
pattern. 

Highpass  filtering  in  the  frequency  domain  is  used  to  accen- 
tuate edge  detail  in  the  spatial  domain  due  to  the  fact  that 
edges  contribute  greatly  to  high-frequency  components.   A 
highpass-filter  transfer  function  is  one  which  attenuates 
frequency  components  within  a  specific  radius  about  the 
origin  defined  by  the  cutoff  point.   The  cutoff  point  is 
chosen  by  examining  the  energy  distribution  in  the  frequency 
domain.   This  distribution  represents  the  amount  of  energy 
contained  within  a  certain  radius  about  the  origin. 

There  are  a  number  of  transfer  functions  which  can  be  used 
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for  highpass  filtering.   The  difference  among  transfer  func- 
tions is  described  by  the  manner  of  transition  between  the 
minimum  and  the  maximum  values.   A  Butterworth  highpass 
transfer  function  was  used  to  filter  the  image  in  Figure 
5.2-4.   This  transfer  function  is  given  by 

1 


H(u,v)   =   ^n" 

1  +  0.414[D0/D(u,v)]zn 


(5.2-4) 


where  DQ  is  the  cutoff  point,  n  is  the  filter  order,  and 
D(u,v)  is  the  distance  from  the  origin  of  the  pixel  at  the 
coordinate  (u,v).   Note  that  the  value  of  the  transfer  func- 
tion at  DQ  is  0.707  rather  than  0.5.   The  energy  distribu- 
tion revealed  that  9.24%  of  the  energy  was  contained  within 
a  radius  of  15  pixels  about  the  origin.   Conseguently,  the 
image  was  filtered  using  a  Butterworth  highpass  filter  of 
order  n=l  with  a  cutoff  point  of  DQ=15. 

After  completion  of  the  filtering  operation,  the  filtered 
image  is  transformed  back  to  the  spatial  domain  using  the 
inverse  DFT.  The  translation  process  is  then  repeated  to 
restore  the  image  to  its  original  phase. 

The  values  for  the  filtered  image  were  found  to  range  from 
-173  to  86  which  is  well  outside  the  acceptable  range  of 
gray  level  values.   The  histogram  for  this  image,  shown  in 
Figure  5.2-5,  indicates  that  a  very  small  amount  of  the  gray 
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levels  are  contained  below  the  value  -80.   As  a  result,  the 
range  was  truncated  to  the  interval  from  -80  to  86  by- 
thresholding  as 


g(x,y)  = 


f(x,y)    if   f(x,y)  >  -80 


-80       otherwise  (5.2-5) 

Translation  to  an  acceptable  range  was  obtained  by  addition 
of  an  offset  constant,  specifically, 

g1(xfy)  =  g(x,y)  +  80.  (5.2-6) 

Figure  5.2-6  shows  a  picture  of  the  resulting  image.   The 
transformed  gray  levels  range  from  0  to  166.   Examination  of 
the  histogram,  shown  in  Figure  5.2-7,  indicates  that  most  of 
the  gray  levels  are  contained  within  the  interval  from  0  to 
140.   The  distribution  generated  a  mean  value  of  81  with  a 
variance  of  385. 

The  gray- level  range  was  expanded  through  contrast  stretch- 
ing.  A  parabolic  transformation  function  was  used  so  that 
larger  values  would  be  affected  to  a  greater  extent.   The 
transformation  function  is  expressed  as 


"         0   if      f(x,y)  <  c1 


g(x,y)  =  < 


255[f(x,y)  -  c1]2   c1  <  f(x,y)  <_  ch 
(ch  -  c,)2 

255  f(x,y)  >  ch 

(5.2-7) 
where  c,  and  c,  are  the  cutoff  points.   Figure  5.2-8  graph- 
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ically  depicts  the  transformation  function. 


e(* ,y) 


2^5 


f(x,y) 


Figure  5.2-8  Parabolic  contrast  stretching 
transformation  function 


The  cutoff  points  c,=0  and  c,=140  were  chosen  for  the  para- 
bolic transformation.   Notable  improvement  in  the  ridges  of 
the  far  left  and  right  portions  of  the  image  is  observed  in 
the  resulting  image  shown  in  Figure  5.2-9.   Examination  of 
the  histogram,  shown  in  Figure  5.2-10,  indicates  that  the 
distribution  is  skewed  towards  the  darker  end  of  the  gray 
scale.   An  increase  in  the  mean  value  to  91  is  a  result  of 
the  expanded  range  of  gray  levels;  however,  the  mean  is 
small  due  to  the  skewness  of  the  distribution.   Contrast 
improvement  is  supported  by  an  increase  in  the  variance  to 
1598. 

Histogram  equalization  was  performed  on  the  image  in  Figure 
5.2-9  to  redistribute  the  gray  levels  uniformly.   An  in- 
crease in  the  degree  of  contrast  is  evident  in  the  equalized 
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image  pictured  in  Figure  5.2-11.   A  substantial  improvement 
in  ridge  detail  is  apparent  in  the  lower  portion  of  the 
image  where  a  majority  of  the  degradation  originally  oc- 
curred.  The  histogram,  shown  in  Figure  5.2-12,  indicates 
that  the  distribution  is  indeed  more  nearly  uniform.   Evi- 
dence of  uniformity  is  reinforced  by  a  mean  value  of  128. 
An  increase  in  the  variance  to  5408  is  representative  of  the 
considerable  contrast  enhancement. 

Smoothing  operations  were  used  to  remove  noise  from  the 
egualized  image  as  well  as  produce  uniform  ridges.   The 
image  in  Figure  5.2-11  was  smoothed  through  neighborhood 
averaging  using  the  3x3  spatial  mask  shown  in  Figure 
3.2-3.   A  picture  of  the  smoothed  image  is  shown  in  Figure 
5.2-13,  with  the  corresponding  histogram  plotted  in  Figure 
5.2-14.   This  operation  resulted  in  a  slight  decrease  in  the 
mean  and  variance  to  127  and  3928,  respectively,  which  is 
typical  of  smoothing  operations.   Nevertheless,  the  degree 
of  contrast  is  guite  acceptable. 

For  the  purpose  of  comparison,  the  egualized  image  in  Figure 
5.2-11  was  also  smoothed  using  a  3  x  3  median  filter.   The 
resulting  image  is  pictured  in  Figure  5.2-15.   To  some 
extent,  the  image  is  less  blurred  than  its  averaged  counter- 
part.  Comparison  of  the  histogram,  shown  in  Figure  5.2-16, 
with  Figure  5.2-12  indicates  a  slight  decrease  in  contrast, 
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which  is  corroborated  by  a  reduction  in  the  variance  to 
4740.   This  reduction  is  not  as  great  as  that  of  the  averag- 
ing operation,  which  indicates  a  greater  degree  of  edge 
preservation.   Essentially,  the  mean  was  unchanged  by  the 
median  filtering  operation. 


Edges  in  the  smoothed  image  of  Figure  5.2-13  were  enhanced 
via  the  gradient  operation.   Detection  of  edges  in  the  hori- 
zontal and  vertical  directions  was  achieved  through  the 
use  of  five-level  masks  denoted  as  G  and  G  ,  respectively, 
and  shown  in  Figure  5.2-17.   The  composite  gradient,  denoted 
as  G,  was  computed  as 


G  =  max[  |Gx|,  |Gy|] 


(5.2-8) 


-1    -2    -1 
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Figure  5.2-17  Spatial  masks  used  to  compute  horizontal 
and  vertical  gradients 


Refinement  of  the  edges  was  achieved  by  subtracting  the 
absolute  value  of  the  Laplacian  gradient  from  the  composite 
gradient  [4].   The  Laplacian  is  a  second-order  derivative 
computed  using  the  spatial  mask  shown  in  Figure  5.2-18. 
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Figure  5.2-18   Spatial  mask  used  to  compute 
Laplacian  gradient 


Gt  = 


Once  the  modified  gradient  is  computed,  an  appropriate 
threshold  is  chosen  so  that  only  relevant  edge  detail 
remains.   Generally,  thresholding  is  performed  using  the 
relationship 

0  if  G  <  t 
255      G  >  t 

(5.2-9) 
The  thresholding  operation  was  modified  for  the  enhancement 

operation  so  that 

r 

1  if  G  <  t 
255       G  >  t 

(5.2-10) 
The  histogram  for  the  modified  gradient,  shown  in  Figure 
5.2-19,  is  helpful  in  selecting  a  proper  threshold.   Values 
for  the  gradient  were  found  to  range  from  0  to  838.   Using 
trial  and  error,  a  threshold  of  300  was  chosen  for  the  gra- 
dient Gt«   Figure  5.2-20  shows  a  picture  of  the  gradient 
after  thresholding. 


Gt  = 
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Edges  were  enhanced  in  the  smoothed  image  of  Figure  5.2-13 
using  the  relation 

f(x,y)    if  G  <  t 


g(x,y)  = 

255         6  >  t 

(5.2-11) 

where  f(x,y)  is  the  smoothed  image  and  t  is  the  threshold 

value  chosen  for  G.  in  Equation  (5.2-10).   This  operation 

was  accomplished  by  multiplying  the  smoothed  image  in  Figure 

5.2-13,  denoted  as  f(x,y),  by  the  threshold  gradient  G  of 

Equation  (5.2-10)  and  truncating  all  values  above  255.   The 

enhanced  image,  shown  in  Figure  5.2-21,  reveals  greater 

contrast  and  significantly  sharper  ridge  detail  as  a  result 

of  the  operation.   The  mean  and  variance  increased  to  139 

and  5136,  respectively. 

A  block  diagram  detailing  Procedure  One  is  presented  in 
Figure  5.2-22. 
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Original  Sample  f(x,y) 


LINEAR  CONTRAST  STRETCH 


SPATIAL  TRANSLATION 


DFT 


i 


HIGHPASS  FILTER 


INVERSE  DFT 


SPATIAL  TRANSLATION 


NORMALIZATION: 
THRESHOLD  AND  OFFSET 


POLYNOMIAL  CONTRAST  STRETCH 


HISTOGRAM  EQUALIZATION 


SMOOTHING 


GRADIENT  THRESHOLDING 


t 

Enhanced  Image  g(x,y) 
Figure  5.2-22  Block  diagram  for  Procedure  One 
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5.3  PROCEDURE  TWO 

The  second  procedure  executes  smoothing  operations  prior  to 
filtering.   Smoothing  the  image  provides  a  degree  of  homoge- 
neity so  that  uniform  ridges  are  presented  to  the  filtering 
operations.   Smoothing  can  be  accomplished  using  either 
neighborhood  averaging  or  median  filtering. 

Neighborhood  averaging  was  used  to  smooth  the  image  in  Fig- 
ure 5.1-1.   This  operation  was  executed  using  the  spatial 
mask  of  Figure  3.2-3.   Noticable  improvement  in  ridge  defi- 
nition as  well  as  a  meaningful  reduction  in  noise  is  evident 
in  the  smoothed  image  shown  in  Figure  5.3-1.   The  histogram, 
shown  in  Figure  5.3-2,  produced  a  slight  decrease  in  the 
mean  value  to  134.   However,  the  variance  actually  increased 
to  837  which  is  unusual  for  a  smoothing  operation. 

Figure  5.3-3  shows  the  result  of  smoothing  the  image  in  Fig- 
ure 5.1-1  using  a  3  x  3  median  filter.   Median  filtering 
produced  a  lesser  degree  of  uniformity  than  neighborhood 
averaging  as  seen  in  the  varying  degree  of  contrast  in  the 
image.   The  accompanying  histogram  is  shown  in  Figure  5.3-4. 
This  distribution  yielded  a  decrease  in  the  variance  to  725 
while  the  mean  was  unchanged  at  135. 

The  smoothing  operation  was  followed  by  highpass  filtering 
in  order  to  enhance  the  ridge  detail.   As  mentioned  previ- 
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ously,  the  smoothed  image  is  translated  prior  to  transforma- 
tion to  the  frequency  domain  so  that  the  origin  is  at  the 
center  of  the  frequency  plane.    Figure  5.3-5  shows  the 
image  of  Figure  5.3-1  in  the  frequency  domain  displayed 
using  the  relation  in  Equation  (5.2-3).   The  effects  of 
smoothing  are  observed  when  this  image  is  compared  to  the 
image  in  Figure  5.2-4.   A  reduction  of  high-frequency  compo- 
nents along  with  the  clustering  of  low-frequency  components 
is  characteristic  of  image  smoothing.   Removal  of  the  spikes 
in  the  second  and  fourth  quadrants  is  evidence  of  a  reduc- 
tion in  noise. 

Filtering  was  performed  using  an  exponential  highpass  fil- 
ter.  The  transfer  function  for  the  exponential  highpass 
filter  is  given  as 

H(u,v)  =  exp[-0.347(D0/D(u,v))2n] 

(3.3-1) 
where  n  is  the  filter  order,  DQ  the  cutoff  frequency,  and 
D(u,v)  the  distance  from  the  origin  of  the  pixel  at  the  spa- 
tial coordinate  (u,v). 

Calculation  of  the  energy  distribution  revealed  that  a 
radius  of  8  pixels  about  the  origin  contains  10.335%  of  the 
energy.   Thus,  the  filter  parameters  were  chosen  as  Dn=8  and 
n=l.   After  filtering,  the  image  was  transformed  back  to  the 
spatial  domain  using  the  inverse  DFT  followed  by  a  spatial 
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translation  to  recover  the  original  phase  information. 

Post-filtering  gray  level  values  ranged  from  -189  to  87. 
Examination  of  the  histogram,  shown  in  Figure  5.3-6,  indi- 
cates that  very  few  pixels  contain  values  less  than  -50. 
Therefore,  the  image  was  truncated  using  a  threshold  of 
t=-50.   An  offset  constant  of  c=50  was  used  to  translate  the 
gray  level  range  to  the  interval  from  0  to  137.   The  normal- 
ized image  is  pictured  in  Figure  5.3-7.   A  plot  of  the 
normalized  histogram,  shown  in  Figure  5.3-8,  produced  a  mean 
of  46  and  a  variance  of  193. 

Histogram  egualization  was  executed  in  an  effort  to  redis- 
tribute the  gray  levels  more  uniformly  throughout  the  range 
from  0  to  255.   A  picture  of  the  equalized  image,  shown  in 
Figure  5.3-9,  exhibits  a  noticable  improvement  in  contrast. 
Considerable  improvement  is  visible  in  reducing  the  amount 
of  degradation  in  the  lower  portion  of  the  image  so  that 
individual  ridges  are  distinguishable.   Most  importantly, 
faint  ridges  in  the  far  left  and  right  portions  of  the  image 
are  appreciably  enhanced.   The  equalized  histogram,  plotted 
in  Figure  5.3-10,  yielded  a  mean  of  130  and  variance  of 
5435. 

A  block  diagram  of  Procedure  Two  is  detailed  in  Figure 
5.3-11. 
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Original  Sample  f(x,y) 


SMOOTHING: 
Neighborhood  Averaging 
Median  Filter 


SPATIAL  TRANSLATION 


HIGHPASS  FILTER  | 


INVERSE  DFT 


SPATIAL  TRANSLATION  | 


NORMALIZE: 
THRESHOLD  AND  OFFSET 


HISTOGRAM  EQUALIZATION 


Enhanced  Image  g(x,y) 


Figure  5.3-11  Block  diagram  for  Procedure  Two 
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Figure  5.  1-1   Sample  recovered  after  two  hour: 
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Figure  5.  2-2   Result  of  linear  contrast  stretching 
applied  to  image  in  Figure  5.  1-1 
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Figure  5.2-4  Result  of  transforming  image  in  Figure  5.2- 
2  to  the  frequency  domain  and  displaying 
using  Equation  (5.2-3) 
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Figure  5.2-6 


Result  of  highpass  filtering  and 
normalization  applied  to  image  in  Figure 
5.2-2. 
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Figure  5.2-9    Result   of    parabolic   contrast   stretching 
applied  to  image  in  Figure  5.2-6 
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Figure  5.2-11   Result  of   histogram  equalization  applied  to 
image  in  Figure  5.2-9 
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Figure  5.2-13   Result  of   neighborhood  averaging  applied  to 
image  in  Figure  5.2-11 
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Figure  5.2-15   Result  of   applying  3   x  3   median  filter  to 
image  in  Figure  5.2-11 
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Figure  5.2-20   Gradient  of   image  in   Figure   5.2-13   after 
thresholding 
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Figure  5.2-21  Result  of  edge  enhancement  through  gradient 
thresholding  applied  to  image  in  Figure  5.2- 
13 
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Figure  5.3-1    Result   of   smoothing   through   neighborhood 
averaging  applied  to  image  in  Figure  5.1-1 
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Figure  5.3-3    Result  of   3  x   3  filtering  applied  to  image 
in  Figure  5.1-1 
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Figure  5.3-5 


Result  of  transforming  image  in  Figure  5.3-1 
to  the  frequency  domain  and  displaying  using 
Equation  (5.2-3) 
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Figure  5.3-7  Result  of  highpass  filtering  and 
normalization  applied  to  image  in  Figure 
5.3-1 
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Figure  5.3-9    Result  of   histogram  equalization  applied  to 
image  in  Figure  5.3-7 
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CHAPTER  6:  SUMMARY  AND  CONCLUSIONS 

The  objective  of  this  project  was  to  enhance  degraded  fin- 
gerprint images  using  digital  image  processing  technigues. 
Improvement  in  image  guality  can  be  utilized  by  criminal  in- 
vestigators to  identify  unrecognizable  fingerprints. 

Prints  for  this  project  were  recovered  from  human  skin  using 
the  modified  iodine-silver  plate  method.   Considering  that 
prints  on  human  skin  deteriorate  after  a  short  period  of 
time,  the  recovery  time  of  latent  prints  placed  on  human 
skin  becomes  critical.   An  increase  in  this  recovery  time  is 
a  conseguence  of  the  enhancement  of  degraded  prints. 

Ridges  are  the  primary  features  to  be  enhanced  in  finger- 
print images.  The  criteria  for  enhancement  of  ridge  struc- 
ture involves  distinguishing  individual  ridges  to  the  point 
where  features,  such  as  ridge  endings  and  bifurcations,  are 
recognizable.  Ridges  must  be  distinct  from  one  another  in 
order  to  precisely  establish  the  relative  positions  of  such 
features. 

Combinations  of  spatial-domain  and  freguency-domain  enhance- 
ment technigues  were  used  to  improve  the  print  images. 
Results  show  that  freguency-domain  technigues  work  well  in 
enhancing  ridge  detail.   Highpass  filtering  was  guite  suc- 
cessful at  bringing  out  ridge  detail  from  areas  which  were 
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either  faint  or  obscured  by  degradation.   Spatial-domain 
techniques  are  best  suited  to  contrast  enhancement.   Histo- 
gram equalization  and  contrast  stretching  were  both  useful 
in  sharpening  the  contrast  between  ridge  and  valley  regions. 
Smoothing  operations  were  utilized  in  order  to  remove  noise 
and  produce  uniformity. 

Two  procedures  were  developed  to  enhance  fingerprint  images 
recovered  through  use  of  a  modification  of  the  iodine-silver 
plate  method.   Procedure  One  combines  contrast  stretching 
and  highpass  filtering  in  order  to  enhance  ridge  detail. 
Histogram  equalization  was  added  to  improve  contrast,  fol- 
lowed by  smoothing  operations  for  additional  uniformity  in 
ridge  areas.   Edges  were  further  enhanced  through  gradient 
thresholding.   Procedure  Two  executes  smoothing  prior  to 
highpass  filtering.   This  is  performed  to  remove  noise  and 
establish  uniformity  in  the  ridge  areas  for  increased  effi- 
ciency in  the  filtering  operation.   Filtering  was  followed 
by  histogram  equalization  so  that  gray  levels  were  redis- 
tributed throughout  the  normalized  range  for  increased  con- 
trast.  Both  procedures  showed  significant  enhancement  of 
ridge  detail.   The  most  appreciable  improvements  occurred  in 
areas  where  ridge  detail  was  either  faint  or  obscured  by 
degradation. 

One  disadvantage  of  the  procedures  is  the  amount  of  trial 
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and  error  involved  in  choosing  suitable  parameters  for  the 
operations.   It  would  prove  beneficial  to  program  the  pro- 
cedures so  that  parameters  are  chosen  automatically.   This 
could  be  accomplished,  for  instance,  by  calculating  the 
cumulative  distribution  function  from  the  probability  dens- 
ity function  to  select  proper  threshold  and  cutoff  values 
for  spatial-domain  technigues.   Compass  gradients  could  also 
be  programmed  to  detect  edges  oriented  in  any  direction  with 
greater  efficiency.   In  the  freguency  domain,  the  energy 
distiribution  provides  insight  into  useful  cutoff  points. 
Automation  of  the  procedures  would  allow  technicians  to  im- 
plement the  enhancement  technigues  in  a  practical  environ- 
ment.  Further  work  along  these  lines  would  provide  a  use- 
ful application  of  this  project. 
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APPENDIX  A:  EQUIPMENT  LIST 

Digital  Equipment  Corporation  VAX  11/750  digital  computer 

Grinnell  GMR  270  Series  High  Performance  Frame  Buffers  and 
Image  Processing  System 

Mitsubishi  Model  #C3922  LPK  High  Resolution  Color  Television 
Monitor 

Image  Resource  Videoprint  5000 

Bencher,  Incorporated  Illuma  System  Quartz  Light  Table 

Hewlett-Packard  7475A  Plotter 

Olympus  35  mm  Camera 

Ilford  ASA  400  35  mm  film 
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APPENDIX  B:   COMPUTER  LISTINGS 

Contained  within  this  appendix  are  computer  listings  of  some 
of  the  enhancement  routines  included  in  the  programs  COLOR 
and  ENHANCE.   The  programs  contain  a  considerable  number  of 
menu-  and  screen-drivers,  which  are  omitted  for  the  sake  of 
brevity.   Only  routines  which  are  relevant  to  the  thesis  are 
listed. 

The  routines  require  certain  header  files  which  contain 
variable  and  array  declarations  for  the  registers  and  work- 
ing arrays.   These  header  files  are  included  in  the  routines 
at  compile  time.   Most  of  the  arrays  are  common-block  arrays 
so  that  each  routine  affects  the  same  memory  locations. 

The  header  files  are  listed  first  followed  by  the  routines. 
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tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  3nd  Computer  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename*,  color. h 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

ROUTINE:        color. h 


DESCRIPTION:    This  is  a  header  file  used  to  declare  the 
common  block  of  d3ta  which  is  used  for 
inteSer*2  i«3de  files.   This  header 
declares  the  two-dimensional  image  array 
and  the  length  of  the  rows  and  columns. 


DOCUMENTATION 

FILES:         None. 


AUTHOR*.         Ronald  A.  Gifford 
DATE  CREATED:    February  17»  1988 

revisions: 

MttMttttttttttttttttttUtttttttttttttttttttttttttttttUtttttttt 


inteser*2  image*  rd>  3r»  bl 

dimension  imaged  :512>  1 1512),    rd(H512f  lt512)i 
+  ir<lJ512fi:!512)i  bl(l!5i2tlt512) 

inteser*4  row*  col»  rrow>  rcolf  3row»  3col» 
+  brow,  bcol 

common  ms3e>  rd»  Sr>  blr  row?  col»  rrow»  rco!»  drow* 
+         Scol»  brow*  bcol 

inteaert4  MAX_PELS>  HAX..LINFS 
parameter  <MAX..PELS=512»  HAX_I_INES=512) 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

* 

Department  of  Electrical  and  Computer  Engineering        * 
Kansas  State  University  * 

* 
VAX  FORTRAN  source  filename:  density. h  * 

t 
tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt* 


routine: 


density >h 


DESCRIPTION:    This  header  defines  the  co»»on  block 
variables  for  the  density  slicing 
routines* 


documentation 
files: 


author: 


Ronald  A.  Bifford 


DATE  CREATED:   February  26*    1988 


revisions: 


ttttttt*ttttttttttttt*ttttttttt*tttttttttttttttttttttttttttttt*tt 


inte3er*2  rlo»  rhi»  blo»  bhit  «Jlo>  3hi>  ylo»  yhi>  »lo» 
+  »hi>  do?  chi»  wlo»  whi 

coiion  /density/  rlo»  rhii  blot  bhi»  dlo»  shi>  ylo>  uhi» 
+  »lo»  ahii  clor  chi>  wloj  whi 
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ttttttttttutttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  and  Co»puter  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename!  hist»h 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttnttttttt* 


routine: 


hist 


DESCRIPTION:    This  is  a  header  file  which  declares  the 
•einry  used  for  the  histogram  parameters. 


documentation 
files: 


author: 


Ronald  A.  Gifford 


DATE  CREATED:    September  26»  1987 


revisions: 


ttttttttttttttt*tttttt*tt*tttt*tti1tttt1/.ttttttttttttttttttttttttttt 


realt4  histogram^  »in>  isxi  x.-axis 
dimension  histodr3m<0:255) >  x...a>:is<0:255) 


common  /hist/  histogram  min>  max*  x_Bxis 
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ttttntttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  and  Computer  EnSineerind 
Kansas  State  University 

VAX  FORTRAN  source  filename?  mask.h 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

ROUTINE:        mask.h 


DESCRIPTION:    This  is  a  header  file  for  the  sP3ti3l 
•ask  routines* 


DOCUMENTATION 

FILES:         None. 


AUTHOR:         Ronald  A.  Gifford 


DATE  CREATED:   December  9»  1987 


revisions: 


*tttttttttt*ttttttttttttt***ttt*tttttttttttttttttt**ttttt*ttttttt 


intesler*4   size 
real*4   mask 

dimension  maskd  J7»l  J7) 
common  /mask/mask •    size 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

t 

Department  of  Electrical  and  Computer  Engineering        * 
Kansas  State  University  * 

* 
VAX  FORTRAN  source  filename!   tath.h  * 

ttttttttttttttttttttttttttttttttttuttttttttttttttttttttttttttttt 


ROUTINE! 


•ath«h 


DESCRIPTIONS    This  is  a  header  file  which  defines  so«e 
co»»onlw  used  aatheaatical  paraneters* 


DOCUMENTATION 

FILES:         None, 


author: 


Ronald  A.  Gifford 


DATE  CREATED!    August  23*  1987 


revisions: 


ttt*tttttttttttttttttttttttttttttttttttttlttttt**t*t*tttttttt*ttt 


real*4  PI 

parameter  (PI  =  3. 141592&535897932) 


real*8  DPI 

parameter  (DPI  =  3.1415926535897932D00) 
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tttttttttnttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Depsrtment  of  Electric3l  and  Computer  Engineering 
Ksnsss  Stste  University 

VAX  FORTRAN  source  filensmeJ  resistor. h 

ttttttttt*ttttttttttttttttttfttttttttttttttttt*t*ttttttt$ttttttt 


routine: 


description: 


register. h 


This  is  3  he3der  file  used  to  declsre  the 
common  block  of  6sta   which  is  used  for 
complex  im3S(?  files.   This  hesder  declsres 
the  two-dimonsional  i»33s  array  and   the 
length  of  the  rows  snd  columns. 


DOCUMENTATION 

FILES:  None. 


author: 


Ronsld  A»  Gifford 


DATE  CREATED:    August  25*  1987 


revisions: 


ttttttttttttttttttttttttttttttttttttttttttttttttttnuttttttttttt 


com?»lex*8  i»3Se»  Y>  stor3ge»  stor3fie2 
dimension  i»3de( I :512»1 :512)»  Y< 1 :512>1 :512) t 

+  stor3Se(i:512fi:512)»  stor3Se2(l J512> 1 :512) 

inteSer*4  rows?  columns?  Yrows>  Ycols?  stolrows>  stolcols» 

+  sto2rowsf  sto2cols 


common  imssei  Y>  storage?  stor3Se2»  rows?  columns?  Yrows? 
+         Ycols?  stolrowsj  stolcols*  sto2rowr>?  sto2cols 

inteSerM  MAX_PELS?  MAX-LINES 
psrsmeter  <MAX_PELS=512?  MAX..LINES=512) 


tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Dersrtaent  of  Electrical  snd  Computer  Engineering 
Ksns3s  State  University 

VAX  FORTRAN  source  filenaae:  Butterworth_his:hP3Ss.for 

ttttutttttttttttttttttttttttttttttttttttttttttttttttttttttttttt* 


routine: 


subroutine  But terworth_hiShP3ss( cutoff? 
order »e»ph3sis) 


DESCRIPTION:    This  routine  filters  3  two-diBension3l 

iasste  3rrsy  using  sn  nth-order  Butterworth 
hiahp3ss  filter. 


DOCUMENTATION 

FILES:         None. 


arguments: 


y 


include: 


return: 


routines 
called: 


author: 


cutoff  resl*4 

Cutoff  freauency. 

order  resl*4 

Order  of  the  filter. 

emphasis        real*4 

High-f reauency  emphasis  constant 


register .h 


lB3Sf?<  ) 


None. 


Filtered  image 


Ronsld  A.  Gifford 


DATE  CREATED:    September  7>  19B7 


revisions: 


**ttttttttttttt*ttt*ttt*tttttttttttttttt*ttttut*ttttttttttttt*tt 
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subroutine  Butterworth_hidhr3ss(cutof f  > order >e«ph3sis) 

******************************** 

Vsriable  declarations   * 

********* t ********************** 

implicit  none 

include  'CGIFFORD. HEADS] redister.h' 
re3l*4  cutoff*  Di  D0>  Hv  e*i>h3sis 
int.eSer*4  i»  J»  order.  N2>  BVfil 

******************************** 

Set  e»rh3sis  factor     * 

******************************** 

if  (flaS  .EQ.  0)  then 

eifhaois  -  1.0 
else 

eaphssis  =  0.0 
endif 

******************************** 

Filter  the  usae        * 
******************************** 

DO  =  sart(2.0*(cutoff**2)> 
N2  -  coluams/2 
srsl  =  N2  -  1 
do  J  =  1»  arsl 

do  i  =  1>  rows 

D  =  sart(1.0*((i-N2)**2  +  (J-N2)**2)> 

H  =  1.0/(1.0  t  0.414*(<D0/n)**(2*order>)) 

H  ■  emphasis  4  H 

inasedjJ)  =  c»r lx(K»0.0)*i»a3e(i » J) 
enddo 
enddo 

do  i  =  1»  srsl 

D  =  sart(1.0*((i-N2)**2)) 

H  ■  1.0/(1.0  +  0.414*((D0/n)**(2*order))) 

H  =  emphasis  +  H 

imase(i»N2)  =  H*i»ss«e(  i  »N2) 
enddo 

iB3Se(N2»N2)  =  eaphasis*in3ac<N2>N2) 
arsl  -  arsl  +  2 
do  i  =  arsl»  rows 

D  =  sart(1.0*((i-N2>**2)> 

H  =  1.0/(1.0  +  0.414*((D0/D)**(2*order)>) 

H  ■  e»phssis  +  H 

i«3Se(iiN2)  =  H*iii3Se(ifN2) 
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enddo 

do  J  =  srfil*  columns 
do  i  =  1»  rows 

D  =  sart<1.0*<(i-N2>**2  +  (J-N2)t*2)) 
H  *  1.0/(1.0  •»•  0.414*<<D0/D)**(?*order)>) 
H  =  esphasi?  +  H 
i»3Se(itJ)  =  H*i»aae<  i  »J) 
enddo 
enddo 

return 
end 


8/: 


tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  and  Computer  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename!  bright-.trans.for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


arguments: 


include: 


return: 


routines 
called: 


author: 


subroutine  bright_t pans ( cutoff > value) 


This  routine  filters  an  iaa3e  file  with 
the  cutoff  point  corresponding  to  a 
brightness  level.   The  resulting  imsSe 
turns  3ll  graw  levels  greater  than  the 
cutoff  to  the  transforation  V3lue. 


DOCUMENTATION 

FILES:         None. 


cutoff  real*4 

The  desired  cutoff  point  for  the 
filter  routine. 

value  realt4 

The  transforation  value. 

register »h 


inage( ) 


None. 


Transformed  image. 


Ronald  A.  Gifford 


DATE  CREATED:    July  7*  1987 


revisions: 


mtmttmmmmtmttttmmmmmummmmmmt 


subroutine  brisht_tranr,<cutof f t v3lue) 
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t 

* 
tttttttttttttttttttttttttttttttttttttttt* 

t  Variable  declarations.         * 

ttttttttttttttttttttttttttttttttttttttttt 

t 

t 

* 

* 
* 
* 


implicit,  none 

include    TGIFFQRn. HEADS] register .h' 

int.eser*4    ii    J 

realt4  cutofft  value 

do  J  =  It  col  nans 
do  i  -  It  rows 

if  (real(i«33e<it J))  ,6E,  cutoff)  then 

iaade(itJ)  -  c»plx< value* 0.0) 
endif 
enddo 
enddo 

return 
end 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

t 

Department  of  Electrical  3fid  Computer  Engineering        * 
Kansas  State  University  t 

* 
VAX  FORTRAN  source  filename?  dark_.trans.for  * 

* 
ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


ARGUMENTS: 


include: 


return: 


routines 
called: 


author: 


subroutine  d3rk_trans(cutof f ) 


This  routine  filters  3n  lmsse  file  with 
the  cutoff  point  corresponding  to  a 
brightness  level.   The  resulting  image 
turns  the  3r3y  levels  below  the  cutoff  to 
the  tr3nsfort3tion  value. 


DOCUMENTATION 

FILES:         None. 


cutoff  real*4 

The  desired  cutoff  point  for  the 
filter  routine. 

value  real*4 

The   transformation  value. 


resiister.h 


ima-Se( ) 


None. 


Transformed  im3_f<?i 


Ronald  A.  Gifford 


DATE  CREATED4.    July  7»  1987 


revisions: 


mtmmmtttmmtmmmmmtmttmmtmmmtsttt 
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subroutine  dark_tr3ns( cutoff »vslue) 
* 

ttttttttttttttttttttttttttttttttttttttttt 
t  Variable  declarations.  * 

*********** ********************** t ******* 
* 

* 
* 
* 
* 


implicit  none 

include    TBIFFORD.HEADSDredister.h' 

inteSerM   i»    J 

real*4  cutoff*  value 

do  J  ■  If  colu»ns 
do  i  =  1)  rows 

if  (real(ifjase(i» J))  .LT.  cutoff)  then 

i»a3e(i*J)  -  c»plx(v3lu{?;0.0) 
endif 
enddo 
enddo 

return 
end 


tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electricsl  and  Computer  Engineering 
K3nsas  St3te  University 

VAX  FORTRAN  source  filensse:   dft2D.for 

***«mmmmm*mtm»mmmm*«mt*«t»t**tt  ***»***»» 


routine: 


description: 


subroutine  dft2D(fl3£) 


This  routine  computes  the  DFT  for  3 
complex  input  seouence.  The  input 
seauence  nust  be  3  two-di*ension3l  seusre 
3rrsy  with  row  3nd  colu«n  lengths  which 
sre  3  rower  of  two. 


DOCUMENTATION 

FILES?         None. 


arguments: 


include: 


return: 


routines 
called: 


author: 


flsS    INTEGER*2 

FIbS  to  deternine  trsnsfor*  direction. 

0  Forwsrd 

1  Inverse 


reSister.h 
•sth. h 


lnssef  ) 


None. 


Trsnsfor*ed  i»sSe. 


Ronsld  A.  Riffard 


DATE  CREATED4.    Austust  17  ►  1987 


revisions: 


«mm«mt*ms**t*mmtm*»t«t«*»***ttmt*m*«mm*mt 


subroutine  dft2D(flsS) 


39 


* 
to****************************** 

*  Variable   declarations        * 

************  <t********t*****<m* 

* 

iiplicit  none 


include  ' TGIFFORD. HEADS! register .h' 
include  'CGIFFORD.HEADS.lmath.h ' 

INTE6ER*2  tlsA 

INTEGER*4  i»  J»  k»  iter>  n»  n2>  n  »_N2»  >tf  z 

REAL*4   sisn.  srUt    nult 


C0MPI.EX*8  W.  TeaPf  work 

dimension  work<15512) 
* 

ttttttttttttttttttttttttttttttttttttttttttttttttt 
t  Calculate  the  nunber  of  iterations.     * 

ttttttttttttttttttttttttttttttttttttttttttttttttt 
* 

iter  ■  0 

n  =  rows/2 

do  while  (n  .GE.  1) 
n  =  n/2 
iter  ■  iter  +  1 

enddo 
* 

**************************************** 

*  Check  for  DFT  direction*        * 

mm*********************************** 
t 

if  <flasi  ,EQ.  0)  then 

si an  =  -1.0 
else 

sisn  =  1.0 
endif 
* 
*****************************:*.***: 

*  Column  transform,      * 

*******************************  ** 
t 

do  x  ■  li  columns 

********************************* 

*  Beain  iteration.        * 

********************************* 

n2  =  rows 
* 

do  i  =  If  iter 
n  =  n2 
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n2  =  n/2 

•ult  =  PI/(1.0*r>2) 
* 

ttttttttttttttttttttttttttttttttntttttttt 

t  Calculate  the  Multiplier*       * 

tmtttsmtsmtmmmmmmttm* 

* 

do  J  =  1»  n2 

3T3  =  (J-l)tault 

W  =  c«Plx(cos(sr3) >si2n*sin(3rs) ) 
* 

ttttttttttttttttttttttttttttttttt 

t  Calculate  partitions.    # 

********************************* 

t 

do  k  =  n>  rows*  n 
■  =  k  -  n  +  J 
»_N2  =  H  n2 

Teiip   =   i»3sle(iD»x)    -   i«3Se<i>_N2>x) 
i«3sfe(»»x)    -    ilft*e(»»x)    +    i»3?.e<i»...N2>x) 
i«3Se(n_N2>x)    =  U*Temp 
enddo 
enddo 
enddo 
* 
************************* 

t  Bit  reverssl.    * 

************************* 

t 

n2  =  rows/2 
J  =  1 

n  =  rows  -  1 
do  i  =  1*  n 

if    <i    fl.T.    J)    then 
Teap    -   i*3<fe(Jrx) 
i«33e(J>x)    -    iasge(i>x) 
iastfedfx)   =  Trap 
endif 
k   =  n2 

do  while  (k  .LT,  J) 
J  ■  J  -  k 
k  =  k/2 
enddo 
J  =  J  +  k 
enddo 
enddo 
* 

************************* 
*  Row  trsnsfora.  * 
************************* 

* 

do  x  =  It    rows 

do  7.    =   1»  coluans 

work  (z)   =   iB3ate(xrz) 


91 


enddo 
t 

ttttttttttttttttttttttttttttttttt 

t  Be3in  iteration*        * 

ttttttttttttttttttttttttttttttttt 
t 

n2  =  coluans 
* 

do  i  ■  1»  iter 

9 


n  -   n 
n2  =  n/2 

■ult  =  PI/(1.0*n2) 
t 

tttttttttttttttttttttttttttttttntttttttt 

t  Calculate  the  aultiplier.       * 

ttttttttttttttttttttttttttttttttttttttttt 
* 

do  J  =  .t  >  n2 

ars  =  (J-l)*ault 

W  =  c»plx(cos(3r2)fsifin*sin(ar£!)) 
* 

ttttttttttttttttttttttttttttttttt 

*       Calculate  partitions.    # 

ttttttttttttttttttttttttttttttttt 


do  k  =  n»  coluansi  n 

■  »  k  -  n  +  J 

»_N2  *   ■  +  n2 

Teap  =  work(n)  - 

work( 

a..N2) 

work(a)  -  work(a) 
work(a_M2)  =  W*Tc 
enddo 

+ 

WO 

rk(a_N2> 

* 

enddo 
enddo 

ttttttttttttttttttttttttt 

t                Bit  reversal.    * 

ttttttttttttttttttttttttt 
t 

n2  =  coluans/2 

J  «  1 

n  =  coluans  -  1 
do  i  =  lj  n 

if  (i  .LT.  J)  then 

Teap  =  work (J) 

work(J)  =  work(i) 

work(i)  *  TeiBF 

endif 

k  =  n2 

do  while  (k  ,LT.  j) 

J  =  j  -  k 
k  =  k/2 
enddo 

J  =  J  +  k 
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enddo 
» 

do  z  =  1 »  colums 

i»3Se(>c»z)  =  work(z) 

enddo 
enddo 
* 

ttttttttttttttttttttttttt 
t  Scslina        * 

ttttttttttttttttttttttttt 
t 

Temp   =  c»pl>:(1.0/(1.0*row5)>0»0) 
do  J  =  1»  columns 

do  i  =  1»  rows 

i»3Se(i?J)  =  i»sse< i > J)*Te»p 

enddo 
enddo 


return 
end 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttutttt 

Department  of  Electrical  and  Computer  Engineering 
K3HS35  St3te  University 

VAX  FORTRAN  source  filena»e:  energy. for 

ttttttttttnttttttttttttttttttututtttttttttttttttttttttttttttt 


routine: 


description: 


documentation 
files: 


arguments: 


include: 


return: 


routines 
called: 


author: 


function  energy< radius) 


This  routine  calculates  the  energy 
distribution  for  3  specific  radius  in  the 
freouency  doaain. 


radius  real*4 

Specified  radius 


register.h 


real*4  Percentage  of  energy  contained 
within  radius. 


None. 


Ronald  A.  Gifford 


DATE  CREATED:   February  28»  19B8 


revisions: 


mmmmtmtmtmmtmmmtmttmmmmtmtmtt 


resl*4  function  energy( radius) 


******************************** 
Variable  declarations   * 

******************************** 


9k 


implicit  none 

include   TGIFFORO.HEADSDreiiister.h' 

inteSerM   i»    J 

real<4  radius*  total*  teir>  rO 


ttttttttttttttttttttttttttttttttttttttttt 

%  Compute  energy  distribution     * 

ttttttttttttttttttttttttttttttttttttttttt 
t 

energy  =  0.0 
total  =  0.0 
do  J  =  1»  colunns 
do  i  =  li  rows 

tesp  =  cabs( iaiticd i J) ) 

tot3l  ■  total  +  te«p 

rO  =  sQrt(real((256-i)tt2  +  <256-J)**2) ) 

if  (rO  ,LE.  r3dius)  thsn 

energy  =  energy  +  temp 
endif 
enddo 
enddo 

energy  =  (energy/total )tl00 .0 
return 
end 
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tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electric3l  and  Computer  Engineering 
Ksnsns  State  University 

VAX  FORTRAN  source  filen3»e:  eou.for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


function  eGuO 


This  routine  nen1rul3t.es  1  512  x  512  i»sge 
file  by  eouslizing  sccordin.g  to  the 
histogram  of  the  origins!  i»33e> 


documentation 
files: 

None  . 

arguments: 

None . 

return: 

integer*2 
0 

1 

iasge< ) 

include: 

register .h 
hist.h 

routines 
called: 

histnom( ) 

author: 

Ronald  A. 

Error  flag 

SUCCESS 

Arrsy  not  normalized 

Enuslized  image, 


DATE  CREATED:    Jsnusry  6i  1988 


revisions: 


ttttttttttttttttttttttttttttttMttttttttttttttttttttttttttttttttt 

integer*2  function  eauO 

m*t  ***************  *********** 
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*       Variable  declaration    * 

tttttttttttitttttttttttttntlttt.it 
t 

implicit  none 


include  TGIFFORD. HEADS: register. h' 
include  'CGIFFORn.HFAPSIJhist.h' 

real*4  ea(0t255)»  px 


inteser*2  err>  histnom 
inteser*4  it   J>  Sr3y 
* 

ttttttttttttttttttttttttttttttttl 
t  Compute  histosrsro       * 

ttttttttttttttttttttttttttttttttt 
1 

err  =  histnorm( ) 
if  (err  »EQ.  1)  then 
eau  =  1 
return 
else 

eau  =  0 
endif 
* 

ttttlttttttlttttttttltltttUttltttllt  1111 
t  Co»pute  uniform  distribution    * 

ttttttttttttttttttttttttttttttttttttttttt 
t 

px  =  0.0 

do  i  =  0»  255 

px  =  px  +  histoSram(  i ) 
ea(i)  =  255*px 
enddo 
* 

Itttttttttttttttttttttttttttttttt 
t  Eaualize  imsae  array    * 

ttttttttttttttttttttttttttttttttt 
t 

do  J  =  1>  columns 
do  i  =  1»  rows 

gray  =  nint(  real  ( ims<de(i  >J) ) ) 
iaaao(i»J)  »  cmpl>:(eu(2ray)  »0.0) 
enddo 
enddo 


return 
end 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttntt 

Department  of  Electrical  and  Computer  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename!  exponential-.hiah.for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


ARGUHENTS: 


include: 


return: 


routines 
called: 


author: 


subroutine  exronenti al_hiah( cutoff > 
order»e*phasis) 


This  routine  filters  a  two-di«ension3l 
i»3de  array  usifMf  3n  nth-order  exponential 
hiahpass  filter. 


DOCUMENTATION 

FILES?  None. 


cutoff  real*4 

Cutoff  freouency. 

order  real<4 

Ord?r  of  the  filter. 

enphasis        real!4 

Hiah-f renuency  eaphasis  constant 


register  .h 


inaae< ) 


None 


Filtered  iiiaae 


Ronald  A.  Gifford 


DATE  CREATED:    September  7»  1987 


revisions: 


t*tttttttt*ttttttttt*t1Httt*tit**tt*t*tt*tttttttt*tttttttt*tttttt 
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subroutine  exponent i ft] _hi ah (cutoff* order > eaphasis) 

******************************** 

Variable  declarations   * 

ttttttttttttttttttttttttttttttt* 

implicit  none 

include  'EGIFFORD.HEAPSDreaister .h' 

real*4  cutoff*  D*  P0»  H*  eaphasis 

inteser*2  flag 

inte3er*4  i>  j?  order*  K'2»  srstl 

******************************** 
Set  eaphasis  factor     t 

******************************** 

if  (flas  .EQ.  0)  then 

eaphasis  =1.0 
else 

eaphasis  =  0»0 
endif 

ttttttttttttttttttttuttttttttn 

Filter  the  iasae       * 

********************** * ********* 

DO  =  sart(2.0*(cutoff**2)) 
N2  =  coluans/2 
ar<3l  =  N2  -  1 
do  J  *  If  ar«!l 

do  i  =  1>  rows 

D  =  sart(1.0*((i-M2)**2  +  <J-N2)**2)) 

H  ■  exp(-0t347*((i:i0/D)**(2*order))) 

H  ■  eaphasis  +  H 

iaase(i»J)  =  caplx(H»0.0)*iBBSe(i>J) 
enddo 
enddo 

do  i  =  1 >  arsl 

D  =  sart(1.0*<(i-N2)**2>> 

H  =  exr<-0.347*((D0/D)**(2*order))) 

H  =  eaphasis  +  H 

i»a3e(i»N2)  ■  H*xm8J!e(i»N2) 

enddo 

i  masse  (N2.N2)  =  e»phasis*in3S:e(N2iN2) 

aral  -  aral  r  2 

do  i  =  arsflj  rows 

D  =  sart(l»0*((i-N2)**2>) 

H  =  exp<-0.347*((B0/D)**(2*order))) 
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H  =  emphasis  +  H 
iaaSe<iiN2)  ■  Hti»S*«(ifN2) 
enddo 

do  J  ■  3ral»  columns 
do  i  =  1 i  rows 

D  =  sart(1.0*<<i~N2)*»2  +  <J-N2)**2) > 
H  ■  exp(-0.347*<(D0/D)**(2*order>)) 
H  =  e»ph3sis  +  H 
i»3Se<i>J)  =   H*i»sse( i » J) 
enddo 
enddo 

return 
end 
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tttttttttttttttttttttttttttttttttttttttttttntttttttttttttttttttt* 

Department  of  Electrical  and  Computer  Engineering 
Kansas  State  University 


VAX  FORTRAN  source  filename?  histabs.for 
tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


subroutine  histabsO 


This  routine  calculates  the  histogram  for 
an  it3de  file.   The  range  of  values  is 
divided  into  256  bins. 


DOCUMENTATION 

FILES:  None. 


arguments: 


include: 


return: 


None. 


register .h 
hist.h 


routines 
called: 


author: 


histogram( ) 
>:_axis( ) 
min  and  max 

None. 

Ronald  A.  Gifford 


Computed  histogram. 
Gray  level  values. 


Minimum  and  maximum  flrsv 
level  values. 


DATE  CREATED:    January  28.  1988 


revisions: 


tt*tttttt*ttttttttttttttttttt**ttttttttttttttttttttttttt*tttttttt 


subroutine  histabsO 
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* 
ttttttttttttttttttttttttttttttttttttttttt 

t  V3ri3ble  declarations  t 

tttttttttttttttttttttttttttttttttttttttt* 
* 

t 

t 


inplicit  none 

include  'CGIFFORD. HEADS] register. h' 
include  ' CGIFFORD. HEADSlhist.h' 

inte3er*4  i>  J»  k 


realty  total >  bin 
total  =  real<colu«ns*rows) 
* 

* 

m*m*tm*tttt**tmt*****m* 

*  Zero  histoarsu  array    * 

mmtmm***************t*»«* 

* 

* 

do  i  =  0»  255 

histodraa(i)  ■  0.0 
enddo 
* 

mm*m**m**«**«m********* 

*  Calculate  bin  width     * 

mmmmmm*************** 

t 

call  ran3e(Bin>aa>:> 
bin  ■  <»3x-»in>/250»0 
x_axis<0)  =  Bin 
do  i  »  it   255 
J  =  i  -  1 

x_sxis(i)  =  x.txis(J)  +  bin 
enddo 
* 

ttttttttttttttttttttttttttttttttt 
t  Coapute  histogram       * 

tttttttttttttttttttttttttttttttt* 
t 

do  J  =  1»  colunns 
do  i  ■  If  rows 
k  =  0 
do  while  (x_axis(k)  .I.T.  real  (i»ase(i>  J) ) ) 

k  =  k  +  1 
enddo 

histossraii(k)    =   histo<3rsni(k)   +   1.0 
enddo 
enddo 
» 

do   i    =  0»    255 

histoara»(i)   =  histoSP3«(i )/total 
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enddo 

return 
end 
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tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electric3l  and  Computer  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename!  histnora.for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


function  hist.nom() 


This  routine  calculates  the  histMrM  for 
an  iiafle  file.  The  iaase  file  is  assumed 
to  he  normalized  between  the  V3lues  0  and 
255. 


DOCUMENTATION 

FILESJ         None. 


arguments: 

None. 

include: 

register. h 
hist.h 

return: 

inteser*2 

0 

1 

histoSraiC ) 

x_sxis( ) 

■in  and  »3X 

routines 
called: 

author: 

Ronald  A.  6 

Error  flas  which 
denotes  acceptable  ranse 
of  values. 

Gr3H  values  within  rans=e 
Bras  values  out  of  r3n3e. 

Computed  histogram. 

Gray  level  values. 

Mini bub  and  bsxjbub 
Sray  leval  values. 


DATE  CREATED4.    January  28»  1988 


revisions: 
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* 
t 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

t 

t 

inteser*2  function  histnormO 

***************************************** 

t  Variable  decl3r3tions  * 

ttttttttttttttttttttttttttttttttttttttttt 

t 

t 

* 


implicit,  none 

include  ' CGIFFORD. HEADS} rc«Jistf?r»h' 
include  'CGIFFORD.HF.ADSIIhist.h' 

integerM  i»  J»  k 


real*4  total 

total  =  l.OKcolumnstrows) 
* 

********************************* 

*  Test  for  normalization  * 

******** ***************** ******** 
* 

call  ranse(min>max) 

if  (Kin  ,LT.  0.0  .OR.  me;-:  .6T.  255.0)  then 
histnorm  -  1 
return 
else 

histnorm  =  0 
endif 
* 
************************* 

*  Scale  x_axis    * 
************************* 

do  i  =  Oi  255 

>:_3>;is<i)  =  reaKi) 
enddo 

* 

********************************* 

*  Calculate  histogram     * 

*********************  ************ 

* 

do  i  =  0>  255 

histodraaU)  =  0.0 

enddo 
* 

do  J  »  it    columns 
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do  i  -  h   rows 

k  ■  nint(  real<  .\«3<3e(i  r  j)  > ) 
histogram (k)  -  histo3ra«(k>  +  1.0 
enddo 
enddo 

do  i  =  0>  255 

histosrsod)  =  histos$ra»(.\  )/tolal 
enddo 

return 
end 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttmtttt 

* 

Depart»ent  of  Electric3l  3nd  Co»Piiter  Engineering        * 

K3ns3s  State  University  * 

VAX  FORTRAN  source  filens»e:  line3r..stretch.for         * 

t 

ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


subroutine  1 ines r_st retch (low» hi gh) 


DESCRIPTION!    This  routine  performs  3  contrast  stretch 
using  3  line3r  tr3nsf or»3tion. 


DOCUMENTATION 

FILES!  None. 


ARGUHFNTS:      low  real*4 

Low  cutoff  point. 


include: 


return: 


routines 
called: 


author: 


high  reaDM 

High  cutoff  point. 


register. h 


i«age( ) 


None, 


Transformed  insgp, 


Ronald  A.  Gifford 


DATE  CREATED:    October  4>  1907 


revisions: 


tttttMttttttttttttttMttttttttttMtMttttttttttttttttttttttttttt 


subroutine  linear_st retch (low  thigh) 


m***m.**mt  ********  ********* 
Vsri3ble  declaration    * 
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ttttttttttt  Ittttt  *************  m 
t 

implicit  none 

* 
* 


include  'CGIFFORD.HFADSJ register, h' 

int.e.<3er*4  i»  J 

resl*4  low>  hiSh>  m  b>  tear 


* 
***************************************** 

*  Calculate  slops  snd  intercept   * 

***************************************** 
* 

m   =  255.0/<hiSh-low) 

b  =  -»*low 
* 
**********  mm***************** 

*  Perform  transforation  * 
********************************* 

* 

do  J  ~  If   colu»ns 
do   i   =   1 t    rows 

if    (resl  (i«isSe(i  r  J) )    ►LTi    low)    then 

ift3*e(irJ)  -  ciF 1m (0.0 5 0.0) 
else   if    (resKinssieU  t  J))    ,GT.    hi<*h)    then 

i»3<Se(i»J)    -  chip  lx(  255.0j  0.0) 
else 

te«p   =  ■*resl (imsde<i t J) )    +   b 
in3£le(i»J)    -  cmpl>:(te«pj0.0) 
endif 
enddo 
enddo 
* 

return 
end 
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tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  and  Computer  Engineering 
Kansas  State?  University 

VAX  FORTRAN  source  filename:  »3sk_.co»pute.for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


subroutine  ■3sk_co«Pute( ) 


This  routine  computes  the  »ask  of  a 
re3l-v3lued  insSe.   The  aask  sire  is 
given  by  the  V3ri3ble  'size'  found  in 
■88k*h«   The  i*a<te  size  is  given  by  the 
variables  'rows*  3nd  ■coluans"  found  in 
re£ister«h. 


DOCUMENTATION 

FILES:         None* 


arguments: 


include: 


return: 


routines 
called: 


author: 


None. 


register. h 
■ask.h 


ifis3e<>  Iiiage  after  «ask   computation. 


c_switch( ) 


Ronald  A.  Gifford 


DATE  CREATED:    December  9t    1987 


revisions: 


tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt* 

subroutine  ■ask_conpute  (  ) 
tt*tttttttttttttttt*tttttttt*t*t 
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*  Variable   declarations        * 

********************************* 

iaplicit  none 


include  '[GIFFORP.HEADSJredister .h' 
include  'CGIFFORD.HF.ADS.Iaask.h' 

intesert4  if  Ji  k»  1»  ik»  Jl»  edse>  lower*  upper 
real*4  sun  teap 


* 
************************* 

*  Coapute  bounds  * 

************************* 
* 

edSe  =  (size-D/2 

lower  =  edse  +  1 

upper  =  rows  -  edSe 
* 
**********************************  *,***.*******»|;*****t*t*t»4**»*** 

*  Coapute  the  task  and  return  value  to  the  iaaainary  part.* 
************************  ********t:****»*t»:************t*.t********* 
* 

do  J  =  lower*  upper 

do  i  =  lower*  upper 
sun  =  0.0 
Jl  =  J  -  edse 
do  1  =  li  size 

ik  =  i  -  edse 
do  k  =  1»  size 

sua  =  sua  +  r<?3Kia3Se(ik>  Jl)  )*aask(k»l ) 
ik  =  ik  \   1 
enddo 

Jl  =  Jl  +  1 
enddo 

ii33e(iu)  =  earlx<  reel  (iaasle<i  >  J) )  >sua) 
enddo 
enddo 
* 
*************************** ************** 

*  Switch  re3l  and  iaastinaru  * 
t%t%%t%t%ttttttttttt%t%t%tt%%%%tt*%tt%ttt 
t 

call  c_switch() 

return 
end 

no 


ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

* 

Department  of  Electrical  and  Computer  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename:  «ean_v3r»for 

ttttttttttttttttttttttttttttttttttttttttttttttttttttttttittttttt 

ROUTINE?        subroutine  nc3n..v3r(«e3n»v3r) 

DESCRIPTION:    This  routine  C3lculst.es  the  aean  3nd 

and  V3risnce  for  3  histogram  of  the  real 
parts  of  3n  iasse  arrsy. 


documentation 
files: 


ARGUMENTS:       «esn     resit* 

Computed  »ean  returned  by  routine. 

var     real>4 

Computed  variance  returned  by  routine. 

INCLUDE:        register. h 
hir.t.h 

RETURN:         None. 


ROUTINES 

CALLED?         None. 

AUTHOR:         Ronald  A.  Gifford 

DATE  CREATED:    December  8>  1987 

revisions: 

*ttttttttttttttttt*t*ttt*tttttttttttttt.tttt*ttttttttttttttttttttt 

subroutine  »ean_.var(iie5nrV3r) 

mm***tm*m************«** 
Variable  declarations    * 
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«****«*t**t**t**mm*mt*tm« 
* 

implicit  none 

* 


include  'CGIFFORD.HEADSlhiEt.h' 
inte=ier*4  i 
resits  team  vsr 


ttttttttttttttttttttttttttttttttttttttttt 
t  Calculate  »ean  and  variance.    * 

mmmmm*t*«***t*tt**t*ttm**m* 

■ean  =  0.0 
do  i  =  0)  255 

■ean  =  dean  +  >;_a>:ic(i  )*histosraii(i ) 
enddo 
* 

var  ■  0.0 

do  i  =  0»  255 

var  =  var  +    ( <>:_3>iis<i  )--iesn)**2)*histosirsm(i ) 
enddo 


return 
end 
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tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  3nd  Computer  Engineering 
Kansas  State  University 

VAX  FORTRAN  source  filename:  median. for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


subroutine  median(size) 


This  routine  performs  median  filtering 
on  two- -dimensional  imagp  arrays. 


DOCUMENTATION 

FILES?  None. 


ARGUMENTS: 


include: 


return: 


routines 
called: 


author: 


sire    integerM 
Mask  size. 


register  .h 


imageO     Filtered  image. 


c_switeh( ) 


Ronald  A.  Gifford 


DATE  CREATED:    February  28*  1988 


revisions: 


mmm*mmmttmmmttmtm*mm*mtmttmttmt 


subroutine  median(size) 

)WM  ******************  ********* 
Variable  declaration    * 

******************************** 

implicit  none 
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* 

include  'CGIFFORn.HEADS.Iredister.h' 
* 

intesfer*4  size*  aed»  lowcr>  upper*    run  cur>>  rl»  rut    rlli 
+  cl»  cu>  i»  Jt   k?  1;  index*  indexl»  index2> 

+  index3 

* 

resl*4  »edc»  »edr»  temp 
dimension  nedc<lt5)>  »edr<U5) 
* 

ttttttttttttttttttttttttt 

t  Initialize      * 

ttttttttttttttttttttttttt 
t 

■ed  =  <size-i)/2  +  1 
rup  =  rows  -  nsed  +  1 
cup  =  colunns  -  «ed  +  1 
* 

ttttttttttttttttttttttttt 
t  Filter  colunns  * 

ttttttttttttttttttttttttt 
t 

do  J  =  «ed»  cup 

do  i  =  medj  rup 

rl  =  i  -  »ed  +  1 
ru  =  i  +  aed  -  1 
el  =  J  -  aed  +  1 
cu  =  J  +  med  -  1 
index  =  1 
do  1  =  cl;  cu 

roedcd)  =  recl(.incfie(rl»l)) 

rll  »  rl  +  1 

index 1  =  1 

do  k  =  rll>  ru 

inde>;2  -  indcxl 

do  while  (real(i«3se(k»U)  »LT. 
+  »edc<index2)  .AND,  index2  ,6E.  1) 

indcx3  -  index2  +  1 
»edc(index3)  =  »edc(  inde>:2) 
index2  ■  index2  -  1 
enddo 

index2  =  index2  +  1 
aedc<index2)  =  r»3l(i»*«e(k» 1) ) 
indexl  =  index!  +  1 
enddo 

■edr(index)  -  nodcfaed) 
index  =  index  +  1 
enddo 
t 

ttttttttttttttttttttttttt 
t  Filter  rows     * 

ttttttttttttttttttttttttt 
t 

index  -  1 
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•edc(l)  =  »edr(l) 
do  1  ■  2>  size 

indexl  ■  index 

do  while  (■edr(l>  ►  LT.  •edc(indexl)  .ANU. 
+  indexl  ,GE.  1) 

inde:;2  *  index  1  +  1 
■edc<index2)  =  »edc( indexl ) 
index 1  =  indexl  -  1 
enddo 

indexl  =  indexl  +  1 
■edc(indexl)  =  nedr(l) 
index  =  index  +  1 
enddo 

i*33e(i»J)  =  c»f1x(  resl(  i«3Se(i>.j) )  »«edc(»ed) ) 
enddo 
enddo 

call  c_switch() 
* 

ttttttttttttttttttttttttttttttttt 
t  Copy  border  values      * 

********************************* 
* 

k  =  nd  -  1 
1  =  rup  +  1 
do  J  =  1 <  k 

do  i  =  It    rows 

tenp  ■  sinas( iBsde(i* J) ) 
i«33e(i>J)  •=  cmplx<tei*p>teap) 
enddo 
enddo 

do  J  =  ied>  cup 
do  i  =  1 i    k 

te«p  =  3im33(i»iD3e(  i » J) ) 
i»3Se(i»J)  =  c»plx( te»p> te*p) 
enddo 
do  i  =  1 i    rows 

te»p  =  3i»s3(i«3£ie(.i  ?J) ) 
i«33e<i>J)  -   c»plx(tf?n>p»te»p} 
enddo 
enddo 

k   =  cup  +   1 
do    J    =    k;    coli.ir.ns 
do  i  »  li    rows 

te»p   =   si»i3S(  i»3.3e(i  r  J)  ) 
i>a3e(i»J)    =   c»plx(te»Pf  terjp) 
enddo 
enddo 
* 

return 
end 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

* 

Depsrt»ent  of  Electric3l  sr\6   Computer  Engineering        * 
Kans3s  State  University  * 

* 
VAX  FORTRAN  source  filename:  rsnge.for  * 

ttttttttttttttttttttttttttttttttttttttttnttttttttttttttttttttttt 


routine: 


description: 


subroutine  ranseCfcinma::) 


This  routine  finds  the  range  of  V3lues 
for  the  real  part  of  3  couple;:  i»3ge 
array. 


DOCUMENTATION 

FILES:         None. 


arguments: 


include: 


return: 


routines 
called: 


author: 


•in     re3l$4 
Returned  minimum  value. 

■  ax      real*4 
Returned  «a>:iiiu»  value. 


register .h 
None. 

None. 

Ronald  A.  Gifford 


DATE  CREATED4.    Septe»ber  24,    1987 


revisions:     12/8/8? 


Changed  mm«um  to  r3nse. 


ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


subroutine  range<«in»i>3>;) 


mm***S*tt***«*************** 
Variable  declarations    * 
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ttttttttttttttttttttttttttttttttt 
t 

iBPlicit  none 

* 

* 
* 


include    '[GIFFORD.HEAIiSJreriister.h' 
resl*4   iim    aaxi    tenp 
intedert2   if    J 


* 

*  Search  for  rande.  * 
m*t***t***m*  ***************** 
* 

Bin  =  res.U  inasteU  »1  )  ) 
■  s>:  ■  »in 
do  j  =  1»  columns 
do  i  ■  1 i  rows 

teftp  ■  real (.\B3de(x > J) ) 
if  (teap  .LT.  Bin)  then 

tin  =  i(?«p 
else 

if  (teap  .  GT.  «?>:)  then 

B3X  =  tc?BP 

endif 
endif 
enddo 
enddo 
return 


end 
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ttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Depsrtient  of  Electrical  and  Computer  Engineering 
Kansas  St3te  University 

VAX  FORTRAN  source  filenaaeJ  slice. for 

tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


description: 


subroutine  sliceO 


This  routine  nerfoms  density  slicing 
using  3  tot3l  of  eight  colors. 


DOCUMENTATION 

FILES?         None. 


arguments: 

None. 

include: 

color,  h 

density. h 

return: 

red() 

blueO 

Sreen( ) 

ROUTINES 

called: 

None. 

author: 

Ronald  A. 

Thresholded  values 
Thresholded  values 
Thresholded  vslues 


DATE  CREATED*.   February  26,  1988 


revisions: 


t***tttttttttttt***tttttttt*tttt*ttttt*ttttttttttttttttttttttttt* 

subroutine  sliceO 
tttttttttttttttttttttt*t**tttt*t 
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*  Variable   declarations        * 

ttttttttttttttttttttttttttttttnt 
t 

implicit  none 


include 

l  T.GIFFORD.HEADSlcolor.h' 

include  'CGIFFORD.HEADSldensita . 

h' 

* 

intedei 

'*4  i»  J 

* 

rrow  = 
rcol  = 
a  row  = 
acol  = 
brow  = 
bcol  = 

row 
col 
row 
col 
row 
col 

* 

do  J  = 
do 

1  >  col 
i  =  1 »  row 
rd(i>J)  ■  0 

3  r  ( i  »  J )  -  0 
bl(i»J)  =  0 

enddo 

enddo 

* 

*«««** ttttttttttttttttttt 

t 

Threshold       * 

* 

tttttttttttttttttttttt 

do  J  = 

1  -  col 

do 

i  =  1»  row 
if  (iaa*r(i»J)  .GE. 

rlo 

.AND. 

+ 

imaasdtj)    .LE. 
rd<i»J)  =  255 

rhi) 

then 

else  if  <  iiri32e(  i  >  J) 

,GE. 

bio 

.AND, 

* 

insstedr  J) 
bl(i»J)  ■  255 

.LE. 

bhi) 

then 

else  if  (i»sSe(i>J) 

.GE. 

.41  o 

.AND. 

+ 

ii»3<?c(  i  *  J) 
Sr(i»J)  =  255 

.LE. 

2hi) 

then 

else  if  (i*3<Je(iiJ) 

.GE. 

wlo 

.AND. 

+ 

illlldi  J) 
Sr(i>J)  =  255 
rd(i»J)  =  255 

.LE. 

yhi) 

then 

else  if  (iB33e(itJ) 

.GE. 

nlo 

.AND. 

+ 

ilticdi  J) 
rd(i»J)  b  255 
bl(i»J)  =  255 

.LE. 

■hi) 

then 

else  if  (ii33e(i»J) 

.GE, 

clo 

.AND. 

+ 

i»33e(i »J) 
Sr(i.J)  =  255 
bl(i»J)  =  255 

.LE. 

chi) 

then 

else  if  <in3s«e(i»j) 

►  GE. 

wlo 

.AND. 

+ 

nasp(i.j) 

.LE. 

whi) 

then 
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Sr(ifJ)  =  255 
rd(itJ)  =  255 
bl<i»J)  =  255 
endif 
enddo 
enddo 
return 
end 
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tttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  3nd  Computer  EnSineerin.3 
Kansas  State  University 

VAX  FORTRAN  source  filename:  sauare_stretch»for 

tttttttttttttnttttttttttttttttttttttttttttttttttttttttttttttttt 


routine: 


subroutine  sau3re_stretch(x0>xl ) 


DESCRIPTION:    This  routine  performs  a  contrast  stretch 
usina  a  simple  souared  polynomial. 


DOCUMENTATION 

FILES:         None, 


arguments: 


include: 


return: 


routines 
called: 


author: 


xO      real*4 

Lower  cutoff  point 

xl      real*4 

Upper  cutoff  point 


register. h 


imaSe( ) 


None. 


Transformed  imase. 


Ronald  A.  Gifford 


DATE  CREATED*.    October  4.  1907 


revisions: 


February  29>  1988   Add  cutoff  points 


mtmmmmttmmmtmtttmttmtmtttmmmttttm 


subroutine  soua re_st retch < xO » xl ) 


tttttttttttttttttttttttttttttttt 
Variable  declaration    * 
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ttttttttttttttttttttttttttttttttt 

t 

implicit  none 
* 

* 

* 


include  'CGIFFORD.HEADS.Iredister  .h' 

inte*er*4  it   j 

reslt4  x0»  xl>  te»p»  const. 


* 

ttttttttttttttttttttttttttttttttt 

t  Perfora  transforation   * 

ttttttttttttttttttttttttttttttttt 
t 

const  =  255.0/<<xl-x0)*t2) 
do  J  =  1>  columns 
do  i  =  1»  rows 

if  ( reel iii23e(i > J)  )  ,LT.  xO)  then 

!■«*«( ir J)  -  cbp1x(0,0»0.0) 
else  if  ( reel ( i»s3e(i > J) )  »GT.  xl)  then 

i»3Se(i»J)  -  carlx(255»0;0,0) 
else 

te«p  =  const*(  (reel  (i«aste(i>  J))-x0>**2) 
i»33e(i>J)  -  c»plx(tei»p»0.0) 
endif 
enddo 
enddo 
* 

return 
end 
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ttttttttttttttttttttttttttutttttttttttttttttttttttttttttttttttt 

Department  of  Electrical  and  Cosruter  Ensineerind 
Kans3s  State  University 

VAX  FORTRAN  source  filename!  translate. for 

**************************************************************** 


routine: 


description: 


subroutine  tr3nslate(> 


This  routine  translates  an  i»3se  file  in 
the  spatial  domain  so  that  the  Fourier 
representation  is  centered  in  the? 
freauencw  sauare.   The  array  mist  be  a 
souare  array. 


DOCUMENTATION 

FILES?  None. 


arguments: 


include: 


return: 


routines 
called: 


author: 


None. 


register. h 


inasie.h 


None. 


Ronald  A*  Gifford 


Translated  image. 


DATE  CREATED:    August  27»  1987 


revisions: 


tttttttttttttttttttttt*t*ttttttttttt**ttttttttt*tttttttttttttMtt 


subroutine  translateO 


******************************** 
Variable  declarations.   * 

*******  ************************* 
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iftplicit  none 

include  'C6IFF0RD. HEADS] register. h' 

inteaerM  if  J 

resltA    rem 


t 

ttttttttttttttttttttttttttttttttt 

*       Translate  i»sae.        * 
ttttttttttttttttttttttttttttttttt 

t 

do  J  =  1»  roluans 
do  i  -   1»  rous 

re«  =  <(i+J>/2.0)  -  <(i  +  .j)/2> 
if  <re»  ,NE.  0.0)  then 

i«3Se(i»J)  -  cbp1>;(-1  .0  »0.  0)*ii33e<  i  u) 
endif 
enddo 
enddo 


return 
end 
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ABSTRACT 

This  thesis  presents  the  results  of  digital  image  enhance- 
ment techniques  applied  to  degraded  fingerprints.   A  brief 
discussion  of  fingerprint  identification  and  recovery  tech- 
niques is  included  for  background. 

Various  image  enhancement  techniques  which  were  employed  are 
explained.  Enhancement  techniques  are  divided  into  spatial- 
domain  and  frequency-domain  operations. 

Two  procedures  which  produced  good  results  are  outlined. 
These  procedures  consist  of  combinations  of  several  enhance- 
ment operations.   Details  of  the  results  obtained  from 
applying  these  procedures  to  a  fingerprint  sample  recovered 
using  the  modified  iodine-silver  plate  method  are  incor- 
porated into  the  discussion. 
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